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MANAGEMENTUITTREKSEL 


titel 

Analysis  &  synthesis  of  (SAR)  waveguide  phased  array  antennas 
(ESA  contract  no.  10134/93/NL/PB) 

auteur(s) 

Ir.  HJ.  Visser 

datum 

februari  1 994 

opdrachtnr 

- 

IWP-nr. 

:  710 

rapportnr. 

:  FEL  -94-B062 

Dit  rapport  beschrijft  methoden  voor  de  analyse  en  synthese  van  golfpijp  phased  array  antennes 
alsmede  enkele  ontwerpen  van  C-band  SAR  dubbel  gepolariseerde  golfpijp  phased  array 
antennes.  Analyse  en  synthese  kunnen  gebruikt  worden  voor  enkel  gepolariseerde 
(rechthoekige)  golfpijp  phased  array  antennes  alsmede  voor  dubbel  gepolariseerde  (vierk ante) 
golfpijp  phased  array  antennes.  Analyse  en  synthese  staan  de  aanwezigheid  van  irissen  in  de 
golfpijp  uitmondingen  en  de  aanwezigheid  van  een  dielektrische  plaat  voor  de  antenne  apertuur 
toe.  Irissen  en  dielektrische  plaat  worden  toegepast  voor  Wide  Angle  Impedance  Matching. 

De  gerapporteerde  werkzaamheden  vloeien  voort  uit  werkpakket  1000  van  ESA  contract  Nr. 
10134/93/NL/PB. 

Begonnen  wordt  met  een  literatuurstudie  gewijd  aan  dubbel  gepolariseerde  golfpijpstralers, 
resulterend  in  de  keuze  van  de  vierkante  open  golfpijp.  Na  een  grondige  beschrijvjng  van  de 
mode  matching  analysemethode  voor  een  oneindige  golfpijp  array  -  eindigheidseffecten 
inbegrepen-,  die  de  basis  vormt  van  al  de  verder  beschreven  analyse-  en  synthese-methoden, 
wordt  de  nauwkeurigheid  van  de  analyse  software  gevalideerd  middels  vergelijk  met  metingen 
aan  twee  gerealiseerde  antennes.  Deze  antennes  bevatten  gecentreerde  irissen  in  de 
golfpijpuitmondingen  en  een  dielektrische  plaat  voor  de  antenne  apertuur  voor  impedantie- 
aanpassingsdoeleinden.  Een  antennesynthesemethode,  gebruik  makend  van  simulated  annealing 
en  downhill  simplex  technieken,  wordt  vervolgens  beschreven  en  verschillende  antenne- 
ontwerpen,  gebaseerd  op  een  enkel  element  in  een  oneindige  array  omgeving,  worden 
gepresenteerd.  Deze  ontwerpen  zijn  gesyntetiseerd  met  en  zonder  de  aanwezigheid  van  een 
aanpassingsstructuur  en  voor  rechthoekige  alsmede  driehoekige  array  roosters.  Vervolgens 
worden  de  ontwerpen  van  subarrays  gepresenteerd  die  aile  voldoen  aan  de  door  ESA  gegeven  en 
in  het  rapport  genoemde  specificaties.  Getoond  wordt  het  grote  belang  van  het  meenemen  van 
de  array  omgeving  in  het  ontwerp  van  een  subarray.  Een  microstrip  patch  golfpijp  exciter  en  een 
microstrip  subarray  voedingsnetwerk  worden  besproken  en  de  diepte  van  de  golfpijpstraler 
wordt  afgeschat.  Gekozen  is  voor  de  realisatie  van  een  rechthoekig  rooster  array  met  golfpijpen 
van  2,5  cm  diepte  zonder  irissen  en  zonder  dielektrische  plaat,  gegroepeerd  in  lineaire  8 
elements  subarrays. 

Deze  realisatie  zal  alleen  dan  plaatsvinden  als  ESA,  na  vergelijk  met  een  concurerende  patch 
antenne  optie  -  aangeboden  door  de  universiteit  van  Karlsruhe  beslist  in  het  voordeel  van 
FEL-TNO.  Deze  beslissing  zal  plaatsvinden  op  review  meeting  1,9  februari  1994. 
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1  INTRODUCTION 


Future  spacebome  SAR  (Synthetic  Aperture  Radar)  systems  should  feature  wider  bandwidth, 
scanning  capabilify  in  azimuth  and  elevation,  and  above  all  they  should  provide  full  polarisation 
information  because  polarimetry  is  one  of  the  most  promising  tools  for  the  interpretation  of 
radar  signatures./ 

Current  advanced  spacebome  SAR  systems  are  not  able  to  acquire  the  full  polarisation  matrix 
and  do  not  have  sufficient  polarisation  purity.  Breadboard  results  from  a  current  state  of  the  art 
system  show  only  17  dB  polarisation  isolation. 

New  radiating  board  concepts  should  therefore  be  investigated  for  providing  higher  polarisation 
purity,  larger  bandwidths  and  low  losses  fulfilling  the  necessary  electrical  and  mechanical 
demands  required  for  space  application. 

This  report  describes  the  development  of  a  high  performance  waveguide  radiating  element  in  an 
array  environment  for  application  in  future  spacebome  active  phased  array  antennas.  Described 
are  the  tasks  as  laid  down  in  workpackage  1 100  due  to  ESA  contract  no.  1 0 1 34/93/NL/PB  ( 1  ]. 
Key  features  of  the  element  in  the  array  environment  are:  Increased  signal  bandwidth,  low 
ohmic  losses,  high  polarisation  purity  and  utilisation  of  low  mass  technology  suitable  for  use  in 
space. 

Specifications  of  the  antenna  are  given  by  [2,  3): 


•  Number  of  ports  per  subarray 

•  Frequency 

•  Bandwidth 

•  Scan  volume 

•  Polarisation 


•  Co  to  cross  polarisation  for  each  scan  direction 
for  each  port 

•  Complex  ratio  of  copolarisations  for  each  scan 
direction  after  calibration  on  boresight  within 

•  Input  reflection  coefficient  "Sjj"  for  each  port 

•  Isolation  "Sj  i+1"  between  2  ports  of  subarray 

•  Maximum  value  of  "active"  reflection 
coefficient  over  scan  volume  for  each  port 

•  Absolute  copolarisation  gain  variations  over 
range  of  polarisation  states 


2  (one  per  polarisation) 
5300  MHz 

100  MHz  (operational) 

±  25  MHz  (margin) 

±  15°  (elevation) 

±  2°  (azimuth) 

Dual  polarisation 
H  and  V  (baseline)  with 
agility  (TBC) 

>  20  dB  (25  dB  goal) 

Modulus  ±  0.3  dB 
Phase  ±  3° 

<  -15  dB 
>30  dB 

<-12.5  dB  (<-  15  dB  goal) 
±  0.4  dB  (±  0.3  dB  goal) 


1  Polarimetry  -  a  technique  used  in  microwave  remote  sensing  -  consists  of  transmitting  horizontally 
and  vertically  polarised  waves,  and  measuring  both  polarisations  of  the  return  signal.  A  target 
polarisation  matrix  is  formed,  and  the  return  from  any  other  incident  polarisation  could  then  be  readily 
synthesized.  This  technique  has  been  used  for  image  enhancement  and  target  discrimination  in 
synthetic  aperture  radars. 
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•  Ohmic  line  loss  in  subarray  :  <  1  dB/m 

•  Aperture  efficiency  ’  :  >  85  %  (>  90  9c  goal) 

In  order  to  fulfill  these  requirements,  the  design  of  such  an  element  in  an  array  environment  has 
to  be  carried  out  by  accurate  modelling  and  not  by  experimental  methods.  Experimental  design 
of  advanced  phased  array  antennas  is  a  very  time  consuming  and  therefore  costly  process  which 
does  not  guarantee  an  optimum  solution.  Especially  when  the  requirements  become  more 
stringent,  acceptable  designs  can  hardly  be  obtained  by  experimental  methods.  Good  theoretical 
solutions  are  needed  to  make  the  antenna  design  simple,  accurate  and  efficient.  Computational 
methods  are  therefore  preferred  because  they  allow  the  optimisation  of  the  various  antenna 
dimensions  and  the  determination  of  tolerances  in  a  simple  way. 

The  emphasis  in  this  report  is  on  the  used  software  tools  for  the  antenna  development  and  the 
related  mathematics  because  of  the  paramount  importance  of  accurate  software  tools  for 
developing  future  spaceborne  SAR  phased  array  antennas. 

This  report  starts  with  a  literature  study  on  waveguide  radiators,  resulting  in  the  choice  of  the 
open  ended  square  waveguide.  After  that,  a  description  of  the  mode  matching  method  used  for 
analysing  a  single  waveguide  element  in  an  infinite  array  environment  will  be  discussed. 

The  reason  for  presenting  this  rather  detailed  description  is  given  by  the  importance  of  this 
method.  All  software  tools  described  in  this  report  rely  heavily  on  this  mode  matching  analysis 
method. 

The  mode  matching  analysis  yields  reflection  and  transmission  coefficients  for  a  single 
waveguide  element  in  an  infinite  array  environment.  Throughout  the  report,  mathematical 
derivations  of  other  antenna  parameters  will  be  given  where  needed. 

After  the  description  of  our  used  analysis  method,  the  report  deals  with  software  validation, 
antenna  synthesis  tools,  first  designs  of  SAR  antennas,  subarray  designs,  the  design  of  a 
waveguide  exciter  and  subarray  feeding  network  and  information  on  the  radiating  plane 
thickness. 


.  .  Embedded  Subarray  Gain  at  Input  Port 

2  Aperture  Efficiency  = - ;-  -c  c — - — r — — - - - 

v  1  4k  Subarray  Aperture 
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2  LITERATURE  STUDY  ON  WAVEGUIDE  RADIATORS 


Of  the  possible  dual  polarised  antenna  elements,  the  microstrip  patch  radiator  and  -  to  a  lesser 
extent  -  the  open-ended  square  and  circular  waveguide  radiator  are  the  most  commonly 
encountered  ones  in  the  literature. 

Other  reported  dual  polarised  elements  are  crossed  dipoles  [4,  5],  printed  log-periodic  antennas 
and  flared  slot  (notch)  antennas  [6j. 

The  problems  associated  with  microstrip  patch  radiators,  i.e.  small  frequency  bandwidth,  high 
losses  and  high  sensitivity  to  tolerances,  both  on  dimensional  and  on  material  parameters,  do  not 
exist  for  waveguide  radiators.  This  is  due  to  the  inherent  wide  bandwidth  characteristics  of  the 
latter.  So  waveguide  radiators  are  strong  candidates  in  the  design  of  a  spacebome  SAR  phased 
array  antenna. 

Waveguide  radiators  known  to  realise  dual  polarisation  are:  Quadruple-ridge  loaded  waveguides 
[7],  Wheeler  lab  square  waveguide  [8],  circular  open-ended  waveguides  [9]  and  square  open- 
ended  waveguides  [10]. 

2.1  Quadruple-ridge  loaded  waveguides 

The  transverse  geometry  of  quadruple-ridged  square  and  circular  waveguides  is  shown  in  figure 
1  18]: 


There  are  two  reasons  why  the  quadruple-ridged  waveguides  are  not  suitable  for  our  purpose. 
First:  Quadruple-ridged  loading  of  square  waveguides  does  not  improve  bandwidth  and  in  the 
case  of  quadruple-ridged  loading  of  circular  waveguides,  the  bandwidth  is  only  moderately 
broadened  [7],  Second:  The  method  of  excitation  of  the  element  is  extremely  difficult  to 
implement  at  high  frequencies,  as  stated  in  [8]. 
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2.2  Wheeler  lab  square  waveguide 

A  sketch  of  the  Wheeler  iab  element  is  shown  in  figure  2  [8]: 


This  Wheeler  lab  element  is  an  unloaded  square  waveguide  with  aperture  dimensions  less  than 
A.j/2,  where  is  the  wavelength  at  the  highest  operating  frequency. 

Due  to  the  reported  poor  cross  polarisation  behaviour  [8]  and  the  mechanical  complexity,  this 
element  is  not  suitable  for  our  purposes. 

2.3  Open-ended  square  waveguide 

Although,  as  stated  in  [9],  circular  waveguides  seem  to  be  the  best  choice,  because  of  reported 
better  matching  characteristics  over  a  scan  volume  of  a  60  degrees  half-angle  cone,  we  prefer 
open-ended  square  waveguides. 

This  because  of  the  relatively  small  scan  volume  required  [1],  because  a  20  percent  bandwidth 
design  has  been  demonstrated  [11],  because  the  ease  in  construction  and  because  of  the  long 
standing  experience  of  TNO-FEL  with  rectangular  waveguides  [12],  combined  with  the 
availability  of  in  house  developed  accurate  software  for  analysing  and  synthesising  open-ended 
square  and  rectangular  waveguide  phased  array  antennas. 
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3  ANTENNA  ANALYSIS  TOOLS 


In  order  to  gain  knowledge  in  the  field  of  polarimetric  phased  array  antennas,  at  TNO-FEL  a 
FORTRAN  code  .for  an  IBM  type  Personal  Computer  is  developed  to  analyse  phased  array 
antennas  consisting  of  rectangular  or  square  cross-section  open-ended  waveguides. 

The  analysis  method  described  in  the  following,  results  in  the  reflection  and  transmission 
coefficients  of  a  single  waveguide  radiator  in  an  infinite  array  environment  as  well  as  the 
transverse  electromagnetic  fields  inside  and  outside  the  antenna  element. 

3.1  Mode  matching 

The  waveguide  phased  array  characteristics  are  determined  by  assuming  that  the  array  is  two 
dimensional  infinite  with  uniform  amplitude  excitation  and  linear  phase  taper,  so  permitting  the 
use  of  a  unit  cell  approach,  which  considerably  simplifies  the  analysis.  The  unit  cell  is  the 
smallest  area  in  which  the  array  can  be  subdivided. 

On  this  unit  cell  a  mode  matching  method  is  applied.  The  field  problem  is  solved  at  the  terminal 
plane  (i.e.  the  waveguide  aperture)  by  applying  the  boundary  conditions  that  the  tangential 
components  of  electric  and  magnetic  fields  are  continuous  across  this  interface. 

First,  the  representations  of  the  tangential  electric  field  when  z  i  0  and  z  t  0  are  matched  to  the 
tangential  electric  field  in  the  aperture  at  z  =  0.  Next  the  representation  of  the  tangential 
magnetic  field  when  z  i  0  is  matched  to  the  one  as  z  T  0,  yielding  an  infinite  system  of  linear 
equations  which  can  be  solved  numerically  by  truncating  it.  This  results  in  the  unit  cell 
reflection  and  transmission  coefficients. 

3.2  Antenna  geometry 

The  array  is  assumed  to  be  two  dimensional  infinite  and  placed  in  the  plane  ?  =  0  of  a  Cartesian 
coordinate  system  x.y,z  (see  figure  3): 


t'igure  3: 


Array  configuration 
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The  array  elements  are  waveguides  of  rectangular  or  square  cross  section.  Centred  irises  (used 
for  impedance  matching)  can  be  placed  in  the  waveguide  apertures.  Because  of  the  inftniteness 
of  the  structure,  edge  effects  are  neglected  and  due  to  the  periodicity  of  the  configuration, 
coupling  effects  are  the  same  for  each  of  the  elements,  so  permitting  the  assumption  of  unit  cells 
into  which  each  element  individually  radiates.  The  problem  has  to  be  solved  for  a  single  unit  cel! 
only.  The  geometry  of  the  unit  cell  is  shown  in  figure  4.  This  unit  cell  contains  the  centre 
element  p=0,q=0  (arbitrary  choice). 


y 

A 


radiating  aperture 
(iris) 

/ 


dielectric  sheet 


t3 

r< 


unit  cell 

waveguide  wall 


Figure  4:  Geometry  of  array  and  unit  cell 

The  grid  is  described  by  the  vectors  s  and  t. 

In  front  of  the  waveguide  apertures,  a  dielectric  sheet  can  be  placed.  Irises  and  dielectric  sheet 
are  used  to  compensate  for  inefficiencies  due  to  coupling  effects.  A  side-view  of  the  unit  cell, 
containing  the  matching  structure,  is  shown  in  figure  5. 
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Figure  5:  Geometry  of  dielectric  matching  structure  and  surroundings 

£,  ,i  =  1,2, 3. 4  denotes  the  relative  electromagnetic  permittivity  in  the  shown  subdomains  of  the 

antenna.  p,;i  =  l,2,3,4  denotes  the  relative  electromagnetic  permeability  in  the  shown 
subdomains  of  the  antenna. 

3.3  Electromagnetic  field  representation 

The  electromagnetic  fields  will  be  described  in  terms  of  modes.  A  complete  knowledge  of  the 
axial  components  of  a  mode  is  sufficient  to  determine  all  field  components  of  this  mode.  As  far 
as  the  field  equations  are  concerned,  the  problem  can  be  separated  in  scalar  problems  in  the 
transverse  plane. 

In  the  remainder  an  exp(-jo)t)  time  dependence  is  assumed  (co=27tf)  and  waves  propagate 
according  to  exp(-jkzz)  for  increasing  z.  kz  is  the  axial  wavenumber. 

For  TE-modes,  the  axial  H-field  components  satisfy  the  two  dimensional  Helmholtz  equation: 


(Vr.Vr  +  *r2)/e=0  (3.3.1) 

For  TM-modes,  the  axial  E-field  components  satisfy  the  two  dimensional  Helmholtz  equation: 
(Vr  •  Vr  +kT2  )e,i  =  0  (3.3.2) 


where  the  subscribt  T  stands  for  the  transversal  part  of  a  parameter.  So  kT  is  the  transversal 
wavenumber. 
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3.3.1  Waveguide  modes 

The  aim  of  this  paragraph  is  to  find  solutions  for  (3.3.1)  and  (3.3.2)  in  the  waveguide  region  of 
the  antenna  and  tQ  derive  expressions  for  the  transverse  fields. 

3.3. 1.1  TE-modes 

As  a  solution  for  (3.3.1)  take: 


hr 


=  / 


d>- 


eosi 


cos 


(3.3.3) 


The  separation  constants  kx  and  ky  are  related  to  the  axial  wavenumber  V ')  by: 

kT:  =*,-+*v2  =C02£0fl0-Y,1|J  (3-3.4) 

with  e0  and  fi()  free  space  electromagnetic  permittivity  and  electromagnetic  permeability, 
respectively. 

I,u  are  mode  amplitude  coefficients, 
imposing  the  boundary  conditions  leads  to: 


+ 


;(m,/t)*(  0,0) 


(3.3.5a) 


(3.3.5b) 


The  axial  wavenumber  is  given  by: 
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(3.3.6a) 


with: 


l 

k0  =  co(e0n0)3 

Further. 


y  <1 ) 
IHmn 


0) 


con( 


■\1H 


1 


y  til 


(3.3.6b) 


(3.3.7) 


3.3.1.2  TM-modes 

Applying  the  same  procedure  as  in  paragraph  3.3. 1 . 1  leads  to. 


k 


i 


(3.3.8a) 


(3.3.8b) 


Further: 


TNO-rspon 


- - - - - - - - -  FEL-94-B062 

Page 

14 


ZEJ"  = 


we 


0 


.  V  (1) 


(3.3.9) 


with  the  axial  wavenumber  given  by  (3.3.6). 


3.3.1.3  Transverse  field  components 

The  transverse  field  components  of  electric  and  magnetic  field  strengths  are  given  by  [12: 


(3.3.10a) 


fy 


(l)! 


(3.3.10b) 


or 


(  .  3 

J 


kr 


±y, 


U) 


(3.3.10c) 


(3.3.1  Od) 


with: 


C((S)s  cos 


(3.3.11a) 
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i,(«- 


sin 


in  E,  +  - 

i  *  2 


(3.3.11b) 


3.3.2  Fields  in  the  terminal  plane 

The  field  distribution  in  the  radiating  aperture  plays  a  vital  role  in  the  analysis  of  the 
configuration  and  should  be  chosen  carefully.  In  van  Schaik  [12],  the  field  is  expressed  as  a 
complete  orthogonal  sequence  of  expansion  functions.  Two  different  expansions  are  described. 
Since  the  expansion  developed  by  Lee  [13]  violates  the  edge  conditions  and  will  limit  the 
permittable  iris  dimensions  (Smolders  [14]),  we  choose  an  expansion  in  Chebyshev  polynomials 
that  satisfy  the  edge  conditions  (but  will  lead  to  more  CPU  time  consumption).  With  this  choice, 
the  tangential  electric  field  in  the  radiating  aperture  of  the  centre  element  is  written  as: 


Ex  ( aperture )  =  £  A*  CTi[^  f~j 


(3.3.12a) 


Ey  ( aperture )  =  £ 


(3.3.12b) 


with: 


C7;(0=[i-(20:p. 


J  cos[iarcsin(2^)];t  =  0,2,4,... 
|-7'sin[tarcsin(2^)];t  =  1,3,5,... 


(3.3.13a) 


STfE)  =  {_Jsin['arcsin(2^)];,‘  =  2,4,6,... 

1  |  cos[i arcsin(24)]t  =  1,3,5,... 

and  A|w,  Bjw  coefficients  to  be  determined. 

3.3.3  Grating  modes  in  the  radiating  domain 


Due  to  the  periodicity  of  the  array  structure,  the  following  relation  holds  for  any  Cartesian  field 
component  in  each  sub-domain  (see  figure  5)  located  in  the  region  z  >  0. 
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{e'M)S(rr  +  r7w),*<n,1(rr  +  rrM)}  =  {«(',)  (rJ.),A,n,"(rr)}exp[>(p4/i+<?T,)]  (3’  '  ' 

'Fj  and  *Ft  denote  incremental  phase  shifts  between  two  adjacent  elements  in  the  directions  ol  s 
and  t,  respectively. 

3.3.3.1  TE-modes 

As  a  solution  of  (3.3.1)  take: 

b;,nf  =  Tn)t  exp (~jkT*rT)  (3.3.15) 


with: 


kT  =  kT  =bt(y¥!k+2pn)  +  b2{yV,+2qn) 


bi  = 


t  xi . 


I,  xs 

b2  =  - - 

*•  r 


(3.3.16a) 


(3.3.16b) 


(3.3.16c) 


S  =  i.  •(sxf)  (3.3. 1 6d) 

Each  set  of  values  p  and  q  specifies  a  wave  component,  called  a  grating  mode,  so: 

V  =7-„“>'exP(-;^.i>)  (3.3.17) 

The  transverse  wave  vector  corresponding  to  the  direction  of  scan  (80,(p0)  is  given  by: 

*r„,  =  *0[sin(eo)cos(<p0)i,  +  sin(eo)sin(<p0)/v]  (3.3.18) 
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The  phase  shifts  can  be  determined  from: 

*¥,  =  kTn  *1  =  k0.s. sin(e0)cos(tp0) 

vfi,  =  k0 .  t.  sin(60  )[cos(Q)cos(cp0 )  +  sin(Q)  sin(cp0 )] 

In  Cartesian  coordinates: 

kTp  =  [(T,  +  2  pit)/  s]ix  +  [('P,  +  2qit)  / 1  sin(Q) -{*¥>  + 2  pn)  /  s  tan(fl)]iv 
Combining  the  results: 

V’,”C’,'*pMv+v)] 

with: 

uno  =  k0  sin(60)cos(<p0)  +  2p7i/  s 

vpq  -  k0  sin(90  )sin(tp0 )  +  2 qiz  / 1  sin(Q)  -  2 pit  /  s  tan(Q) 

The  axial  wavenumber  yw<n)  is  given  by: 


\_ 

iko\  -  “Po2  -  vw2)2  ’ko >  u„o2  +  vw2 

-y(-^o2er„  +«Po2  +  vw2)'  \k02  <  up02  +  vw2 
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(3.3.19a) 


(3.3.19b) 


(3.3.20) 


(3.3.21) 


(3.3.22a) 


(3.3.22b) 


(3.3.23) 


Further: 
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in)  _  . 


ypq  /  wh0;z„  = 


(n) 


3.3.3.2  TM-modes 

Analogous  to  the  procedure  in  paragraph  3.3.3. 1: 

C",!  =  C"’  exp[-X^„o^  +  vMy)] 

Further. 


(3.3.24) 


(3.3.25) 


=  Y  ‘-*/(oefe0;y£i 


(n) 


(3.3.26) 


with  the  axial  wavenumber  given  by  (3.3.23) 

3.3.3.3  Transverse  field  components 

The  transverse  components  of  the  electric  and  magnetic  field  strengths  are  given  by  [10]: 


(nr  _  /_ 


<°Po  V„T„ 


+y 


u  T 

up0l  pq 


(3.3.27a) 


(n)1 


=  {«#**  V"'  +  r„<%  C  )*„  / 

'kv/  ♦‘■w.C  )*„  'V 

-  (* )*„  I K‘ 


(3.3.27b) 


(3.3.27c) 


(3.3.27d) 
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with: 


^sexpj-j(“(^  +  %>’)] 


(3.3.28) 


3.4  Transfer  properties  &  matching 

For  analysing  the  electromagnetic  properties  of  the  dielectric  matching  structure,  in  each 
subdomain,  wave  matrices  and  transfer  matrices  are  introduced.  A  wave  matrix  is  a  column 
matrix  that  contains  the  components  of  the  mode  in  the  direction  of  increasing  z,  and  the 
components  of  the  mode  in  the  direction  of  decreasing  z: 


4/'"r  (z) 

T<«>  (,) 


n  =  2,3,4 


(3.4.1) 


3.4.1  Transfer  matrices 

The  subdomains  and  their  quantities  are  shown  in  figure  6. 
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li  m[  'P  ‘ 2 » ( 2 )]  =  [  77?' 2  *  (  t  2 )]  1  im[  'F ' 2  ’  ( z )]  (3.4.2) 

The  transfer  matrix  of  the  section  t2  <  z  <  T2+T3  ([TR(,,(T3)J)  follows  from: 

lim[^,,,(:)j  =  [r/?°,(t,)]  lim  [4',,,(z)]  (3.4.3) 

The  transfer  matrix  of  the  transition  from  subdomain  2  to  subdomain  3  ([TR2_»3^  f°"ows 
from: 


(3.4.4) 


The  transfer  matrix  of  the  transition  from  subdomain  3  to  subdomain  4  ([TR3 _ >.4])  follows 

from. 


lim  J'F,3'(z)]  =  [r/?M4]:Hmj4'‘4,(2)] 


(3.4.5) 


The  aim  of  the  remainder  is  to  express  |tw<2’  ,7^'  }  in  terms  of  jr^'41  , 7"^  }. 


3.4. 1.1  TE-modes 

(vj/m)]  js  jaicen  to  be: 


* 
H , 


(n)* 

(«)" 


(3.4.6) 


With  the  above,  we  get: 
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lim 

;l0 


H. 

Hz 


(2) 


<2r 


—  f  TRH  1  lim 
1 


H . 

-I 

H. 

L  '« 


<  ■*>' 


(4  >' 


(3.4.7) 


with: 


[  77?W2-.4  ]  =  [7/?i2,(t2  )\TRHi^\T#”{x3)\TRHi-«] 


(3.4.8) 


The  transfer  matrices  [TR<2,(T2)]  and  ITR(J){t)]  are  given  by: 


[TR[n){  xj] 


exp(yYw,n,tJ 

0 


0 


exp(-/irw<n,xj 


;n  =  2,3 


(3.4.9) 


The  transfer  matrices  (TRH2 _ >3]  and  [TRH3_>4]  follow  from  the  continuity  of  the  tangential 

components  of  the  electric  and  magnetic  field  strengths  across  the  interfaces  z  =  T2  and  z  =  T2  + 
T3- 


(3.4.10a) 


lim  {£/”,///”}  =  lim  {et(4,,Ht'4>} 


(3.4.10b) 


This  leads  to: 


[r/?"2_3]=— 1 


2  Y<2) 


yh  <2)  -  yH 


(*) 


k„  <2)  -  ,3) 

Hn  hn 

J7)  +  V,(5) 


(3.4.11a) 


[r/?V4]  = 


2Y<» 


Y„  iy)  +  YH  <4)  Yh  w-Yh  ,4)' 

"no  "bo  Mna  "no 


Yh  m~yH 


(4) 


(4) 


(3.4.11b) 
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Further: 


H2  <2>-  " 

exp 

-y 

(W„0Jr  +  VM>') 

Hz  af 

V2r  exP 

-> 

Mpo*+v„y) 

(3.4.12a) 


lim 

.-it.+T, 


H, 


(4)* 


;4,‘  exp [-AufiX  +  v^y  +  yS' 
0 


(T:+T,))] 


(3.4.12b) 


Now,  the  mode  amplitude  coefficients  {r  <:’  ]  can  be  expressed  in  terms  of  mode  amplitude 
coefficients  }: 

T  t2l‘  =  F  T  (4r  (3.4.13a) 

m  m  n 

T  ,2y  -  G  T  (3.4.13b) 

M  MM 


with: 


Fp*~’ 


( V”  +  V’,K“  +  +  Y„'”l,)]  + 

(V1  -  v‘')=4'h~'!'’’ -y„'”»,)] 


(3.4.14a) 


<V,“ 


(yh„‘2)  ~  YH„ai)(YH„°'  +  -  Yw(,>t3)]  + 

( V2’ + v3,X  v*  -  v^M-yIy,;^: +yw,5,x,)] 


exp[-JfYw,4)(tJ+t3)]/4y'HJJ'yHJ3' 


(3.4.14b) 
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3.4.1.2  TM-modes 

Analogous  to  the  procedure  in  the  previous  paragraph. 


=  F  T 

p<t  r« 


(3.4.15a) 


=  GnJr<l 


(41* 


(3.4.15b) 


with: 


Fr,= 


+v")(  v +zf;4,)^pp(Y(;:,T: +y„o,x,)]+ 


exp[-yyMu,(T:+T))]/4ZiJ:Zi 

(3.4.16a) 


(z£;:-+zf;’’)(zf;',-z£;<')exp[-j(yM<^T:+Yw'’'xJ)] 


exp[-yYw“"(t:  + 1,)]  /  4 Z,  ,:'ZC 
(3.4.16b) 


3.4.2  Matching 

The  representations  of  the  tangential  electric  and  magnetic  field  as  z  X  0  will  be  matched  to  the 
ones  in  the  aperture  at  z  =  0;  the  representations  of  the  tangential  electric  and  magnetic  field  as  z 
1 0  will  be  matched  to  the  ones  in  the  aperture  at  z  =  0.  This  results  in  the  following  equations 
[12,15]: 


(3.4.17a) 
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to/a 


(3.4.1 7b) 


0)^1,, 


?)'■ 


r1"’  -r(" ) 


FiFJkrJ  =  '5LauF\  >F.2-,H\n*0 


(3.4.17c) 


+/J"  )  FtFJkTJ='£BuFi}F„:Lm*0 


(3.4. 17d) 


-C"  )  +  <*, 


mn  mn 


nn  V— «ir  +j-»> 


mn  /  r  H 


Z[-V„,!,>»(F„-Cw)V+<oe(,v„(F(,  +  C„)T„'4,'k|../*,_! 


(3.4. 17e) 


^oupoiFn  +  ^w)]^n^i8 


(3.4.170 


nv 

Fj  to  Fj g  are  defined  in  appendix  A. 

3.4.3  Matrix  representation 

The  matching  of  the  representation  of  the  tangential  electric  field  as  z  i  0  to  the  one  in  the 
aperture  at  z  =  0  yields  a  system  of  the  form: 


T  =  C0)A 


(3.4.18) 


T  being  a  column  matrix  containing  the  free  space  transmission  coefficients. 
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The  matching  of  the  representation  of  the  tangential  electric  field  as  z  t  0  to  the  one  in  the 
aperture  at  z  =  0  yields  a  system  of  the  form: 

R  =  D"'A  +  D'2'l  (3.4.19) 

R  being  a  column  matrix  containing  the  wa%'eguide  reflection  coefficients.  1  is  a  column  matrix 
containing  the  waveguide  excitation  coefficients. 

The  matching  of  the  representation  of  the  tangential  magnetic  field  as  z  i  0  to  the  one  in  the 
aperture  at  z  =  0  yields  a  system  of  the  form: 


U'hT  =  U'2)J  +  U‘'R  (3.4.20) 

In  the  above.  A,  I,  R  and  T  are  column  matrices  formed  by 

B,  }  {/  1~"'  h{/  ,/  h  I.IV  ,A'\T  41  }.  respectively.  The  other  matrices 

are  matrices  of  known  coefficients. 

Substitution  of  (3.4. 18)  and  (3.4.19)  in  (3.4.20)  results  in: 


[vr"-U':,]A  =  fi"’  (3.4.21) 

with: 

W'"  =  Uil)C"  (3.4.22) 

W12’  =  U0'D'2'  (3.4.23) 

B0)  =  t/l4)/  (3.4.24) 

f/<4>  =  t/<2)  (3.4.25) 


The  matching  procedure  leads  to  an  infinite  system  of  linear  equations  with  { A|w,  Bjw)  as 
unknowns.  This  system  can  be  solved  numerically  by  truncating  it.  Once  equation  (3.4.21)  is 
solved  for  A,  the  unknown  column  matrices  T  and  R  are  obtained  upon  substituting  A  in  (3.4.18) 
and  (3.4.19),  respectively. 
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3.4.3. 1  Derivation  of  C(1) 

After  some  mathematical  manipulation  of  equations  (3.4.1 7a, b)  [13],  we  get: 

T  (4)'  —  V1  c  <la,A  4 .V  r  <**> a 

1P1  ~  Z^'-'pqlw  I*  y'pql*  I* 

/h'  Iw 


=2,c„*'m\.+'£cp^wbIw 

/»•  /m 

with: 

c«.'“  *  -(‘r„!  '  +  V)(f„  +  G  Jf  v„F,F,;w  *  0 

<v*' = (v /  f;F,mti0)[(«,l0! + + «„)]"  vyw  *  0 

/  fiKY„n’)[(.v!  +  v)(f^  -  *  o 

Cw-',w  =  I  p.ppi  „'”)[(“«  +  CX^-ep]'1  •v/A:'  *  0 

Or: 

T  =  Ca>A 


FEL-94-B062 


(3.4.26a) 


(3.4.26b) 


(3.4.27a) 


(3.4.27b) 


(3.4.27c) 


(3.4. 27d) 


(3.4.28) 


with: 
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T 4)'-' 


pmax.-r/  max 


l  p  max.  (/max 

*  “  =(4)* 

I  -pmax.-i/max 


(3.4.29a) 


t(4) 

l  p  max, i/max 


;  first  half:  vr  *  0,  second  half:  /  *  0 


(3.4.29b) 


-  p max q  max. 0.0 


jt ./ nux. h  max  C",’l-/>mfix.-4fmax.<l.n  •"  C'^-pmax.-f/max.lmax...  max 


...  C<l‘,'(imax.f(max.(I.O  •"  C(l‘%max.amax.lmax.»  max  C  pmax.^  max.0.0 

C*  ’  =  —  t !«)  731*1 

C  -pmax.-</max.0.0  C  -pmax.-qmax./max.H  max  v*  -  p  max. -y  max. 0.0 


Wl») 

v  pmax.</max./max.n  max 
C-  -pnux.-qmax./max.K  max 


I  - 1  1  Ci )  ^<lu)  *  t 

c  /fmax.ymax.O.O  C  pmax.ffmax.fmax.a  max  t-  pmax.f/max.0.0 

first  quadrant  :  w^O 

second  quadrant  :  1*0 

third  quadrant  :  w  *  0 

fourth  quadrant  :  1*0 


7;d*) 

C-  p  max  ,q  max. /max.*  max 


(3.4.29c) 
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3.4.3.2  Derivation  of  D(1)  and  D(2> 


After  some  mathematical  manipulation  of  equations  (3.4.17c,d)  [13],  we  get. 


or 


=  ^  D,>b)mntwB,w  l mn‘ 

/»•  /** 


(3.4.30a) 


with: 


/u  /H 


(3.4.30b) 


D"u\n,»  =  -j{kTJ 


!(?)'*(" 


^4 


/  nK\ 
b  ) 


FuFn;w,n*  0  (3.4.31a) 


D  mn ih,  — 


mn  V  (  nn  )  }  _ 

-  +  T  |  *5*6 

ay  V  b  J  \ 


— 0  (3.4.31b) 

a  / 


D'-\,,=-j(kr_ <nj") 


JY  rnttV  f  nrc  VlfF  ( 
|V  a  J  V  fo  J  j  3  4  ' 


nm 


FuFn\w, «*0  (3.4.31c) 


(—  k^4;Z,m^0  (3.4.31d) 
V  b  J 


Or: 


R  =  Da)  A +  Di2)J 


(3.4.32) 


with: 
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'00 


tor 


wmax.nmax 


7',r 
/  00 


.(U 

/  mnwx.flnux 


;  first  half:  (m.n)  *  (0,0),  second  half:  m  *  0,«  *  0  (3.4.33a) 


/  = 


00 


1  mmax.rtn 

l  00 


.()> 

I  m  max. n  max 


;  first  half:  (m.n)  *(0,0),  second  half:  m*0,n*0  (3.4.33b) 


rxdt/i  riH«) 

D  0,0 .0,0  *  * '  L)  0.0./ max. h- max 

LJ  o.o.o.o 

rxO/>) 

U  0.0.  /  max.  m  max 

D"'  = 

rx(la) 

L)  m  max  .n  max  .0.0  '**  LJ  m  max.  n  max./  max.  w  max 

—On)  Tt""’ 

U  O.O.O.O  LJ  0.0 ,/ max.  h  max 

LJ  m  max  .n  max  .0.0 

n(l*) 

Lf  o.o.o.o 

r\UM 

•  •  LJ  mmax.rtmax./max.H  max 

*■  U  0.0./  max.  w  max 

—do) 

U  minax.nmax.0.0  Lf  ninu.inuijmai.wiim 

n06’ 

LJ  mmax.nmax.0.0 

T>m 

"  U  mmax.fl  max. /max.  Umax 

first  quadrant  :  w  *  0,«  *  0 
second  quadrant  :  /  *  O.m  *  0 
third  quadrant  :  w*0,m*0,n*0 
fourth  quadrant  :  /  *  0,  m  *  0,  n  *  0 

(3.4.33c) 
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Da)  = 


"-1 

...  0 

0 

O' 

0 

...  -1 

0  ••• 

0 

0 

...  0 

1  ••• 

0 

0 

...  0 

0  ••• 

1_ 

(3.4.33d) 


3.4.3.3  Derivation  of  U(U,  U<2)  and  U(3) 

After  some  mathematical  manipulation  of  equations  (3.4. 1 7e,f)  [13],  we  get: 

Y  U0a)n^TM,A)'  +YJj°a)  +Zt/U  + 

pq  W  m"  mn 

/  <!»'  ,  V77,3u)  ,  7 

>  1/  mnh  +  2-1  U  *  1 


(3.4.34a) 


Tt/"v;r  +Str  + 

w  n  mn  mn 

'ZU'*'-*'  Ij'y  +'£U(ih)mnlJj'y 


(3.4.34b) 


with: 


V''-'«.*<i*oV„(F;+G^)FKF,JkrJ-,l*0 


(3.4.35a) 


(3.4.35b) 
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Ui2u) ^  =  /Y».(,>(~)^3^4  /  kTJ ;/  *  0,m  *  0 


U"u'mr.,.  =  7'toeJ  —  %Fu/kT  2;/*0,m*0 
b 


U*u)mnu=-Uaa)m„u\l*0,rn*0 

77  (3«)  .  '(2u)  .  _ 

U  rnnl w  =  U  mnlw  \  /  *  0,  W  *  0 


^"V  =  -  GJ^F]& /  *  2;  vv  *  0 


t/"V  =coe0M,0(Fw  +  gJ/-7/-8/*, 2;w-*0 


r/(2fc)  _ 

U  mnlw  —  jy„ 


—  l^F;,/*-  2;w*0,«*0 
b  ‘ 


(3.4.35c) 

(3.4. 35d) 

(3.4. 35e) 

(3.4.350 

(3.4. 35g) 

(3.4.35b) 

(3.4. 35i) 


77<2ft) 

U  mn/vt  —  — 7(0£0 


f  W7I 

V  a 


FuFn!kT  2;w*0,n*0 


(3.4.35J) 


1/<3W»*,  =  -U{2h)mniw-,W  *  0,n  *  0  (3.4.35k) 

U°h)mnlw  =  U^mnlw',™*  0,rt  *  0  (3.4.351) 

Or: 

U(')T  =  Ua)I  +  Uti)R  (3.4.36) 
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with: 


tf(lfl) 

U  -/>m*x  ,-q  max.0.0 


I /(la) 

t/  -/imax.-qmax./max.wmax 

77"“' 

lv  -/>max.-qmax.().(» 


77““' 

(_/  -/imax.-^mix./max.H  max 


first  quadrant  :  /  *  0 
second  quadrant :  /  *■  0 
third  quadrant  :  w  *  0 
fourth  quadrant  :  *  0 


U 


(In)  »/(!/>) 

ax.q  max. 0.0  U  -/>  max. -9  max.0.0 


W‘to'  p  max.q  max ./  max .  h  max  u"“  *-/>max.-qmax./max.w  max 

—tin)  JjOM 

U  />  max. q  max.0.0  V  -/>max.-q  max.0.0 


77<1«>  JTito) 

U  pnux.q  max./max.  *max  U  -/»max.-qmax./max.u  max 


i/UM 

U  p  max.v  max.0.0 


r/(l*) 

U  />ma  max. /max. w  max 

77“*' 

ly  p  max.  4  max.0.0 


77"*» 

(7  />max.qmax./max.»*  max 


(3.4.37a) 


(/(:>  = 


r;(2u) 

U  o.o.o.o 

F  f(2«) 

Cy  0.0./ max. u  max 
(2/>) 


uu 


o.o.o.o 


(J{  m max. n max.0.0  U  "  O.O.O.O 


u 


(2«) 


m  max .n  max./  max .  w  max 


y-A2u) 

U  0.0./  max.  w  max 

771 2b) 


m  max  .a  max  ,0.0  U  O.O.O.O 


>  .(2a) 

U  m  max.  n  max.  0.0 


ni2d) 

U  mmax.nmax./max.H  max 
tt(2M 

U  m  max  .n  max.0.0 


|_C/  0,0. /max. h  max 

first  quadrant  :  /  *  0,  m  *  0 
second  quadrant :  /  *  0,  m  *  0,n  *  0 
third  quadrant  :  w  *  0,  n  *  0 
fourth  quadrant  :  w  *  0,  m  *  0,  n  *  0 

(3.4.37b) 


<2*>) 


m  max  .n  max ,/  max  .*  max 


n,2h) 

U  0.0. /max. w  max 


j-ab) 

U  nimax.r  max./max.w  max 
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U 


(21  _ 


|  .(2 a) 

—U  u.o.o.o 

/  /<2</) 

v-7  0.0./  max.  »■  max 

» ji2h) 

—U  «.(),().() 

~~U  0.0./ max.  m- max 


-u{- 


-Uau> 


m  max.  n  max 


Tr  ' 

.0.0  U  o.o.o.o 


m  max. n  max./  max.*  max 

/  #(2M 

(-7  m  max.n  max.  0.0 


//(2M  //“■"' 

—  O'  m  max.n  max./max. h  max  O  0.0,/ max. » 


U  0.0./ max. h  max 

Tj{2b) 

U  o.o.o.o 


mrz(2h) 


first  quadrant  :  /  *  0,  m  /  0 
second  quadrant :  l*0,m*Q,n*0 
third  quadrant  :  u  #  0,  /i  *  0 
fourth  quadrant  :  h-  #  0,  m  *  0,  n  *  0 


I 

O  m  max,  n  max.  0.0 


J-AZu) 

O  m  max. «  maxi  max.  h  max 

Tab, 

U  mmax.nmax.0.0 


C/  m  max.n  maxi  max. m  max 


(3.4.37c) 


Now,  with  equation  (3.4.28)  the  free  space  transmission  coefficients  can  be  calculated  and  with 
equation  (3.4.32)  the  waveguide  reflection  coefficients  can  be  calculated.  The  waveguide 
excitation  coefficients  are  fed  into  column  matrix  I. 

3.5  Finiteness  effects 

Very  large  phased  array  antennas,  i.e.  antennas  consisting  of  thousands  of  elements,  can  be 
analysed  perfectly  well  by  an  infinite  array  theory.  For  medium  sized  phased  array  antennas  (a 
couple  of  hundreds  of  elements),  an  infinite  array  approach  will  do,  but  for  more  detailed 
information,  finiteness  effects  have  to  be  included  in  the  array  behaviour  calculations. 

It  is  because  of  the  above  reason  and  because  of  our  subarray  design  strategy  -  to  be  discussed  in 
chapter  7  -  that  we  have  developed  an  approximate  method  for  calculating  the  behaviour  of  an 
element  in  a  medium  sized  phased  array  antenna. 

Since  direct  calculation  of  finite  array  behaviour  is  very  complicated  and  very  time  consuming 
for  arrays  consisting  of  more  than  say  10  elements,  it  is  expedient  to  look  for  an  approximate 
method. 


3.5.1  Approximate  method 

Use  is  made  of  the  fact  that  reflection  and  coupling  coefficients  are  related  by  Fourier 
transforms  [16].  Reflection  coefficients  as  function  of  scanangles  00  and  (p0  are  calculated  for  an 
element  in  an  infinite  array  environment. 

By  integrating  these  reflection  coefficients  over  the  incremental  phase  shifts  between  adjacent 
elements  in  x  and  y  directions,  coupling  coefficients  between  the  elements  are  approximated. 
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S(myn) 


Jf(  d0^e^e^d(yxdx)d(V^) 


(3.5.1a) 


with: 


V,  =  A0.sin(i30)cos(tp0) 


(3.5.1b) 


Vv  =  A0.sin(t30)sin(<pu) 


(3.5.1c) 


A'0=y  (3.5.  Id) 

The  indices  m  and  n  denote  the  position  of  the  element  within  the  array.  (m,n)  =  (0.0)  is  the 
centre  element. 

Note  -  with  reference  to  figure  4  -  that  for  a  rectangular  grid  array,  the  parameters  (s.t)  and 
(dx,dy)  may  be  interchanged. 

With  a  finite  set  of  so  derived  coupling  coefficients  (the  number  of  coefficients  equal  to  the 
number  of  elements  in  the  finite  array),  reflection  coefficients  are  recalculated,  yielding  the 
approximate  reflection  coefficient  behaviour  of  an  element  in  a  finite  array  environment: 


m  /v 

r(tf0,<p0)  =  ]L  ^S(m,n)e 


-jn\yydx 


(3.5.2) 


with  M  and  N  defining  the  (finite)  size  of  the  array. 

It  appears  that  the  error  made  in  this  approximation  is  much  less  than  the  error  in  applying  an 
infinite  array  approach  as  will  be  shown  in  chapter  4. 


3.5.2  Finite  array  properties 

In  order  to  discuss  the  properties  of  finite  arrays  we  start  with  a  fictive  infinite  array.  The 
dimensions,  with  reference  to  figure  4,  are  given  below,  (these  dimensions  are  taken  from  the 
design  to  be  presented  in  paragraph  6.2). 

•  s  =  t  =  36.26  mm 

•  a  =  b  =  32.94  mm 
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•  c  =  d  =  22.14  mm 

•  T2  =  ° 

•  T3  =  0.77  mm 

•  er  =  6.36 


Coupling  coefficients  relative  to  the  centre  element  (m=0,n=0)  are  calculated  for  finite  arrays 
with  -M<m<M,  -N<n<N  and  M=N.  For  M=31,  the  amplitude  of  the  coupling  coefficients  on 
the  main  axes  of  the  array  are  shown  in  figure  7. 


Figure  7:  Amplitude  coupling  coefficients  relative  to  centre  element,  f -5. 3GHz  TEI0- 

mode  polarisation 

We  see  a  sharp  decrease  in  mutual  coupling  when  moving  away  from  the  centre  element  as 
expected.  Moving  away  more  than  31  elements  from  the  centre  element  does  not  seem  to  be 
meaningful,  since  the  coupling  coefficients  can  not  be  calculated  anymore  with  enough 
accuracy. 

For  M=31,  reflection  coefficients  in  E-,D-  and  H-plane  are  recalculated  with  the  obtained 
coupling  coefficients.  The  results,  together  with  the  infinite  array  reflection  coefficients,  are 
shown  in  figures  8,  9  and  10. 
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Figure  8:  Amplitude  infinite  and  approximate  E-plane  finite  array  reflection,  f=5.3GHz, 

TE  l0-mode  polarisation 


Figure  9:  Amplitude  infinite  and  approximate  D-plane  finite  array  reflection,  f=5.3GHz, 

TE)0-mode  polarisation 
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H-P  LANE 


theta  (degrees) 

-+-  Infinite  *  Units 


Figure  JO:  Amplitude  infinite  and  approximate  H-plane  finite  array  reflection,  f~5.3GHz, 

TE l0-mode  polarisation 

The  above  figures  show  a  near  perfect  match  between  infinite  array  reflection  coefficients  and 
approximate  finite  array  reflection  coefficients  for  M=31.  This  indicates  that  for  this  particular 
array  a  63x63  element  array  can  be  regarded  as  being  infinite. 

In  the  far  scan  region  we  see  a  fluctuation  on  the  finite  array  reflection  coefficient.  This 
phenomenon  is  due  to  the  finiteness  of  the  array  and  is  well  described  in  [16,  17].This 
phenomenon  is  also  observed  in  measurements  on  a  finite  array  [14],  When  more  elements  are 
added  to  a  finite  array,  the  frequency  of  the  fluctuation  increases  and  the  amplitude  decreases. 
To  show  this,  in  figure  1 1  finite  array  reflection  coeffcients  for  different  values  of  M  are 
compared  with  the  infinite  array  reflection  coeffcients  in  the  H-plane. 
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H-PLANE 


0  3 


theta  (degrees) 

-Infinite  *  mmax — 3  *  mm»x  —  7  ^  mmax-IS 


Figure  //:  Amplitude  reflection  coefficient  finiteness  effects.  f^S  JGHz.  TEl0-mode 

polarisation  ( mmax-M ) 
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4  SOFTWARE  VALIDATION 


Our  waveguide  phased  array  antenna  analysis  software  (program  PHASAR)  has  been  tested  by 
calculating  the  unit  cell  reflection  coefficients  for  two  realised  waveguide  phased  array 
antennas.  The  first  antenna  is  CAISSA:  a  planar,  space-fed,  phased  atTay  antenr.a  of  about  850 
radiating  elements  and  operating  in  the  C-band.  The  second  antenna  is  an  EXPAR 
(Experimental  Phased  ARray)  test  antenna,  a  255  element  array  antenna  operating  in  the  X- 
band. 

For  both  antennas  the  dimensions  are  such  that  within  the  scan  range  0-90  degrees,  only  a  single 
propagating  mode  is  radiated  in  free  space. 

4.1  CAISSA  antenna 

The  CAISSA  antenna  is  shown  in  figure  12  [12].  The  dimensions,  with  reference  to  figure  4.  are 
stated  in  table  I . 


Figure  12: 


CAISSA  antenna  without  dielectric  matching  sheet 


The  reflection  coefficients  calculated  with  PHASAR  are  compared  with  the  results  obtained 
from  the  no  longer  available  program  TSCHEBY  [12]  and  results  obtained  from  a  computer 
program  written  by  K.K.  Chan  (Chan  Technologies  Inc.,  Quebec  Canada)  based  on  a  generalised 
scattering  matrix  method  in  combination  with  mode  matching. 

The  calculated  reflection  coefficient  amplitudes  obtained  with  the  three  different  programmes 
are  shown  in  table  2a.  The  related  phases  are  shown  in  table  2b. 


Table  2a: 


Amplitudes  reflection  coefficient  TE]Q-mode  CAISSA  at  5.65GHz 


tp  (degrees) 


0  (degrees)  1  IH-PHASAR  |  1TI-TSCHEBY  |  IH-CHAN 


0 

0 

0 

20 

0 

40 

0 

60 

45 

20 

45 

40 

45 

60 

90 

20 

90 

40 

90 

60 

Table  2b: 


Phases  reflection  coefficient  TEjQ-mode  CAISSA  at  5.65GHz  (degr.) 


0  (degrees) 


Z  T-PHASAR  Z  T-TSCHEBY  Z  T-CHAN 


0 

0 

105.74 

105.50 

103.2 

0 

20 

90.56 

90.25 

92.4 

0 

40 

40.20 

39.90 

55.5 

0 

60 

329.48 

329.60 

327.0 

45 

20 

102.63 

102.21 

100.1 

45 

40 

58.64 

56.79 

73.7 

45 

60 

279.75 

279.87 

289.0 

90 

20 

115.78 

115.31 

108.6 

90 

40 

159.27 

158.96 

127.1 

90 

60 

194.18 

194.71 

134.9 
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From  the  above  it  can  be  concluded  that,  in  comparison  with  similar  purpose  programs,  the 
waveguide  phased  array  antenna  analysis  program  PHASAR  is  working  well.  The  ultimate 
validation  will  of  course  be  found  in  comparison  with  measurements.  This  will  be  discussed  in 
the  next  paragraph. 

4.2  EXPAR  test  antenna 

The  EXPAR  test  antenna  is  shown  in  figure  13  [14],  The  dimensions  are  classified. 


Figure  13:  EXPAR  test  antenna  with  iris  sheet  and  dielectric  matching  sheet  on  top 

4.2.1  Infinite  array  calculations 

For  a  certain  frequency  fg,  the  EXPAR  unit  cell  reflection  coefficients  as  function  of  scan  angle 
8  are  measured  for  H-,  D-  and  E-plane.  These  measurements,  together  with  the  calculated 
reflection  coefficients  are  shown  in  figure  14.  Coupling  coefficients  are  measured  with  a  vector 
network  analyser  and  these  are  fed  into  equation  (3.5.2)  to  get  the  'measured'  reflection 
coefficients. 
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Figure  14:  Amplitude  and  phase  ofEXPAR  unit  cell  reflection  coefficient 
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The  above  figure  is  self-evident,  for  a  basic  antenna  design  the  infinite  waveguide  phased  an-ay 
analysis  program  PHASAR  will  do. 

Next  to  investigate  are  finiteness  effects  for  EXPAR. 

4.2.2  Finite  array  calculations 

Since  we  are  dealing  with  a  wide-angle  waveguide  phased  array  antenna  with  phaseshifters 
applied  on  element  to  element  basis,  finiteness  effects  can  be  analysed  by  integrating  infinite 
array  reflection  coefficients  over  visible  space  only,  as  dicussed  in  paragraph  3.5.1. 

Equation  (3.5.1)  is  now  expressed  as  an  integration  over  i30  and  <p0.  A  phase  correction  factor  is 
added  to  deal  with  the  triangular  lattice.  Furthermore,  the  integration  is  simplified  by  making 
use  of  symmetry  properties  in  the  integration  over  tp0.  Equation  (3.5.1)  is  replaced  by  [18]: 


with: 


,  d,d,  }  } 

S(m,n)  =  4.k0-.—r]  |  r(d0,cp0)cos 


4  n 


o  o 


VdA  m  +  ^ 

*  *  i  2 


cos[  cos  f)0  sin  •d0d$0dy0 

(4.2.1a) 


|  1  n  odd 
P  [O  n  even 


(4.2.1b) 


Equation  (3.5.2)  remains  unchanged. 

Figure  15  shows  H-,  D-  and  E-plane  results  for  a  frequency  fj 
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Figure  IS:  Amplitude  EX  PAR  infinite  and  approximate  finite  array  reflection  coefficient 

The  above  figure  clearly  shows  that  by  approximating  a  finite  array  we  get  very  close  to  reality. 

It  is  thoroughly  demonstrated  now  that  we  have  the  tools  necessary  to  design,  in  a  reliable  way, 
a  SAR  waveguide  phased  array  antenna. 


i 
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5  ANTENNA  SYNTHESIS  TOOLS 


At  TNO-FEL  a  computer  code  is  developed  for  analysing  infinite  arrays  of  open  ended 
waveguides  [12,  45].  The  antenna  geometry  is  shown  in  figure  16.  The  array  lattice  can  be 
triangular  or  rectangular.  The  waveguides  can  be  rectangular  (single  polarisation)  or  square 
(dual  polarisation).  For  impedance  matching  purposes,  thin  irises  can  be  placed  in  the 
waveguide  apertures  and  a  dielectric  sheet  can  be  placed  in  front  of  and  parallel  to  the 
waveguide  apertures. 


radiating  aperture 


dielectric  sheet 


\ 


iris 


A 


waveguide  wall 


z  =  o 


Figure  16:  Open-ended  waveguide  phased  array  antenna  geometry 

Design  parameters  are: 


a 

b 

c 

d 

dx 

dy 

n 

T2 

*3 

Er 


waveguide  width 
waveguide  length 
iris  width 
iris  length 

element  spacing  in  x-direction 

element  spacing  in  y-direction 

triangular  lattice  angle  (90°  for  rectangular  lattice) 

width  air  gap 

width  dielectric  sheet 

dielectric  constant  sheet 
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Given  a  set  of  antenna  design  specifications,  an  antenna  can  be  designed  by  varying  the  design 
parameters  and  using  the  mentioned  analysis  code.  Due  to  the  number  of  design  parameters  and 
the  existence  of  many  local  optima  in  the  solution  space,  this  is  a  very  time  consuming  task  and 
the  finding  of  an  optimum  antenna  design  can  not  be  guaranteed. 

An  optimisation  program  -using  the  analysis  code-  has  been  developed,  that  allows  for 
automatically  generating  an  antenna  design  given  a  set  of  design  specifications  as  input. 

The  optimisation  method  used  is  an  intermixed  combination  of  simulated  annealing  and  the 
downhill  simplex  method  due  to  Nelder  and  Mead  [19],  In  this  chapter,  first  the  basics  of 
simulated  annealing  and  the  downhill  simplex  method  will  be  explained,  after  which  the 
intermixed  combination  of  both  methods  for  continuous  optimisation  will  be  discussed. 

5.1  Simulated  annealing 

Simulated  annealing  methods  are  known  to  have  solved  some  problems  previously  thought  to  be 
practically  insoluble  (for  example  the  traveling  salesman  problem  [20])  and  have  a  direct 
bearing  on  the  problem  of  finding  global  extrema  in  the  presence  of  large  numbers  of 
(undesired)  local  extrema  [21-24],  The  latter  makes  this  method  very  attractive  for  our  synthesis 
problem. 

At  the  heart  of  the  method  of  simulated  annealing  is  an  analogy  with  thermodynamics  [20], 
specifically  with  the  way  that  liquids  freeze  and  crystallise,  or  metals  cool  and  anneal.  At  high 
temperatures,  the  molecules  of  a  liquid  move  freely  with  respect  to  one  another.  If  the  liquid  is 
cooled  slowly,  thermal  mobility  is  lost.  The  atoms  are  often  able  to  line  themselves  up  and  form 
a  pure  crystal  that  is  completely  ordered  over  a  distance  up  to  billions  of  times  the  size  of  the 
individual  atom  in  all  directions.  This  crystal  is  the  state  of  minimum  energy  for  this  system.  If, 
however,  a  liquid  metal  is  cooled  quickly,  it  does  not  reach  this  state  but  rather  ends  up  in  a 
polycrystalline  or  amorphous  state  having  somewhat  higher  energy. 

The  so-called  Boltzmann  probability  distribution: 


Prob(£)  =  exp(-£  !  kT ) 


(5.1.1) 


expresses  the  idea  that  a  system  in  thermal  equilibrium  at  temperature  T  has  its  energy 
probalistically  distributed  among  all  different  energy  states  E.  Even  at  low  temperature,  there  is 
a  chance  of  a  system  being  in  a  high  energy  state.  Therefore,  there  is  a  corresponding  chance 
for  the  system  to  get  out  of  a  local  energy  minimum  in  favor  of  finding  a  better,  more  global, 
one.  The  quantity  k  is  Boltzmann’s  constant  that  relates  temperature  to  energy. 

In  other  words,  the  system  sometimes  goes  uphill  as  well  as  downhill;  but  the  lower  the 
temperature,  the  less  likely  is  any  significant  uphill  excursion. 

In  order  to  use  the  method  described  above  for  optimisation  purposes,  one  needs: 

•  A  description  of  possible  system  configurations 

•  A  random  generator  to  probe  the  solution  field 

•  A  function  E  to  be  minimised  (optimised) 

•  A  control  parameter  T  and  an  annealing  schedule  to  lower  it 
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In  case  of  the  optimisation  of  our  waveguide  antenna,  the  possible  system  configurations  are 
given  in  a  set  of  design  parameter  ranges.  The  function  to  be  minimised  is  the  unit  cell  reflection 
coefficient.  The  control  parameter  T  is  also  the  unit  cell  reflection  coefficient.  Starting  value  of 
T  and  annealing  schedule  are  problem  dependent  and  require  some  experimentation  120).  This  is 
the  major  drawback  of  this  method.  But,  if  for  a  specific  problem  the  annealing  schedule  is 
figured  out,  an  optimum  is  found  that  is  better  than  found  by  any  other  method. 

A  downhill  move  is  allways  accepted;  the  chance  of  accepting  an  uphill  move  is  prescribed  by 
the  Boltzmann  probability  distribution  if  no  additional  information  of  the  minimisation  function 
is  known.  The  optimisation  can  be  speeded  up  by  taking  a  modified  probability  distribution  [21]. 
This  technique  is  known  as  Fast  Simulated  Annealing. 

5.2  Downhill  simplex  method 

The  downhill  simplex  method  is  due  to  Nelder  and  Mead  [19],  A  big  advantage  is  that  this 
method  requires  only  function  evaluations,  not  derivatives.  It  is  an  advantage  because  we  can 
calculate  the  unit  cell  reflection  coefficient  [12,  15],  but  we  do  not  have  an  implicit  function 
description  describing  the  transformation  from  input  parameters  to  reflection  coefficient. 

A  simplex  is  the  geometrical  figure  consisting,  in  N  dimensions,  of  N+l  points  (or  vertices)  and 
all  their  interconnecting  line  segments,  polygonal  faces,  etc. [20],  In  two  dimensions,  a  simplex 
is  a  triangle.  In  three  dimensions  it  is  a  tetrahedron.  Only  nondegenerate  simplexes  -  i.e. 
enclosing  a  finite  inner  N-dimensional  vector  space  -  are  taken  into  account. 

The  downhill  simplex  method  must  be  started  with  N+l  points,  defining  an  initial  simplex.  The 
downhill  simplex  method  now  takes  a  series  of  steps,  most  steps  just  moving  the  point  of  the 
simplex  where  the  function  is  largest  through  the  opposite  face  of  the  simplex  to  a  lower  point. 
These  steps  are  called  reflections,  and  they  are  constructed  to  conserve  the  volume  of  the 
simplex  (hence  maintain  its  nondegeneracy).  When  it  can  do  so,  the  method  expands  the  simplex 
in  one  or  another  direction  to  take  larger  steps.  When  it  reaches  a  "valley  floor",  the  method 
contracts  itself  in  the  transverse  direction  and  tries  to  ooze  down  the  valey.  Note  that  this 
method,  in  contradiction  to  simulated  annealing,  does  not  pretend  to  find  a  global  minimum;  it 
gets  stuck  in  a  local  minimum. 

Basic  moves  for  a  general  3-dimensional  minimisation  problem  are  shown  in  figure  17  [20] 
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Figure  17:  Possible  outcomes  for  a  step  in  the  downhill  simplex  method 


5.3  Continuous  optimisation  by  simulated  annealing 

The  simulated  annealing  method  described  in  paragraph  5.1  is  initially  meant  for  optimising 
multidimensional  discrete  functions.  It  can  be  used  for  continuous  optimisation  but  tends  to 
become  inefficient  as  mentioned  in  [20]  and  as  experienced  by  the  author.  The  inefficiency  is 
due  to  two  phenomena  [20]:  Inefficiency  of  the  random  changes  generator  when  proposing  an 
uphill  move  even  when  local  downhill  moves  exist  and,  most  of  all,  inefficiency  in  convergence 
when  an  optimum  is  approached. 

Therefore,  a  method  proposed  in  [20]  is  used  that  is  an  intermixed  combination  of  simulated 
annealing  and  downhill  simplex.  The  optimisation  starts  as  a  simulated  annealing  method  and 
gradually  changes  to  a  downhill  simplex  method  as  the  temperature  is  lowered.  Thus  at  least  an 
efficient  endgame  is  guaranteed;  when  the  global  optimum  area  is  found  by  simulated  annealing, 
the  downhill  simplex  method  rapidly  isolates  the  optimum. 
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The  system  state  is  described  by  a  simplex  of  N+l  points,  N  being  the  dimension  of  the 
optimisation  problem.  Optimisation  moves  are  as  described  in  paragraph  5.2:  reflections, 
expansions  and  contractions  of  the  simplex.  The  simulated  annealing  is  implemented  in  the 
following  way:  added  is  a  positive,  logarithmically  distributed  random  variable,  proportional  to 
the  temperature,  to  the  stored  function  value  associated  with  every  vertex  of  the  simplex,  and  a 
similar  random  variable  is  subtracted  from  the  function  value  of  every  new  point  that  is  tried  as 
a  replacement  point  The  method  allways  accepts  a  true  downhill  step,  but  sometimes  accepts  an 
uphill  one.  In  the  limit  T  — »  0,  this  algorithm  reduces  exactly  to  the  downhill  simplex  method 
and  converges  to  a  local  minimum. 
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6  BASIC  ANTENNA  DESIGNS  OF  CANDIDATE  SOLUTIONS 


With  the  tools  described  in  the  previous  chapters,  designs  are  generated  based  on  a  single 
element  in  an  infinite  array  environment.  These  designs  will  serve  as  basis  for  the  design  of  a 
subarray.  The  design  of  the  subarray  is  described  in  chapter  7.  The  designs  are  obtained  by 
optimisation  with  and  without  matching  structure  (irises  in  combination  with  a  dielectric  sheet). 
The  optimisation  is  performed  for  rectangular  as  well  as  triangular  array  lattices. 

Optimisation  parameters  are  obtained  from  the  specifications  [2,  3]: 


Number  of  ports  per  subarray 

Frequency 

Bandwidth 

Scan  volume 

Polarisation 

Co  to  cross  polarisation  for  each 
scan  direction  for  each  port 
Input  reflection  coefficient  for 
each  port 

Maximum  value  of  active 
reflection  coefficient  over  scan 
volume  for  each  port 
Ohmic  losses  in  subarray 
Minimum  element  spacing 


2  (one  per  polarisation) 

5300  MHz 

100  MHz  (operational) 

±  25  MHz  (margin) 

±  15  degees  (elevation) 

±2  degrees  (azimuth) 

Dual  polarisation  H  and  V  (baseline) 
with  agility  (TBC) 

>30  dB 

<-15  dB  (<-20  dB  goal) 

<-12.5  dB  (<-15  dB  goal) 

<1  dB/m 
0.5  mm 


Before  optimisation  is  performed,  we  need  to  know  the  number  of  elements  in  one  subarray, 
because  this  parameter  will  determine  the  allowable  element  spacings. 

After  the  size  of  the  subarray  is  determined,  4  designs  will  be  presented. 


6.1  Number  of  elements  in  subarray 

We  start  with  expressing  the  element  distance  as  fraction  of  k: 


del=-  (6.1.1) 

i 

with  i  to  be  determined. 

The  subarray  distances  for  elevation  scanning  and  azimuth  scanning  are  now  given  by: 


dnh  =  M  — ;  elevation 
i 


(6.1.2a) 
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with  M  the  number  of  elements  in  the  subarray  in  vertical  direction. 


dsub  =  N~;  azimuth 
J 

wi*h  N  the  number  of  elements  in  the  subarray  in  horizontal  direction. 


(6.1.2b) 


The  positions  of  the  first  grating  lobe  in  elevation  and  azimuth,  respectively,  are  given  by  1 1 4 j : 


sin(dw  )  =  sin(dri  ) — —  , elevation 

''  "  M 


sin(f)  .  )  =  sin(o  ,  )~  —  \azimuth 

\  u.,  /  \  u.„  /  w 


(6. 1  3a) 


(6.1.3b) 


In  the  above  equations  \>w  and  flu.  are  the  positions  of  the  first  grating  lobes  in  elevation  and 

azimuth,  respectively.  and  are  the  main  beam  positions  in  elevation  and  azimuth, 
respectively. 

The  assurance  that  the  first  grating  lobes  are  positioned  outside  the  scan  volume  for  all  main 
beam  positions  within  this  volume  is  not  enough  for  determining  the  allowable  number  of 
elements  in  the  subarray.  We  have  to  be  sure  that  the  level  of  the  first  grating  lobe  is  sufficiently 
surpressed  within  the  scan  volume.  Therefore  we  look  at  the  angle  where  the  subarray  radiation 
pattern  is  20  dB  down.  With  use  of  [25],  we  find  in  first  order  approximation  for  this  angle: 


0  ,  =  52.03 — (degrees);e/eva/ton 


(6.1.4a) 


d  =  52.03— (degrees),  azimuth 

U‘-!IW  B  \r  c 


(6.1.4b) 


In  table  3,  grating  lobe  angles  (degrees)  and  -20dB  angles  (degrees)  are  given  as  function  of  the 
number  of  elements  in  the  subarray  and  as  function  of  j  (elevation)  or  j  (azimuth)  for  the  main 
beam  scanned  to  the  maximum  angle  ( 1 5  degrees  for  elevation;  2  degrees  for  azimuth): 
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Table  3:  Array  angles  as  function  of  element  spacing  and  number  of  elements 


M(el.)  or  N(az.)  i  (M)or  j  (N)  GL~angle  el.  [  GL-angle  az.  -20dB-angle  ] 


2 

10/5 

10/6 

10/7 

10/8 

48 

35 

27 

21 

75 

53 

43 

36 

52 

43 

37 

33 

3 

■  SH|| 

24 

39 

34 

17 

31 

29 

10/7 

13 

26 

25 

10/8 

9.1 

22 

22 

4 

14 

28 

26 

10/6 

9.1 

22 

22 

10/7 

5.6 

19 

19 

10/8 

3.1 

16 

16 

5 

8  1 

21 

21 

4.3 

17 

17 

10/7 

1.5 

15 

15 

10/8 

- 

12 

13 

6 

4.3 

17 

17 

1.1 

14 

14 

10/7 

- 

12 

12 

10/8 

- 

10 

11 

7 

|  |HH| 

1.5 

15 

15 

1 

12 

12 

10/7 

9.7 

11 

10/8 

8.3 

9.3 

8 

12 

13 

1 

10 

11 

10/7 

8.3 

9.3 

10/8 

7.0 

8.1 

9 

11 

12 

10/6 

8.6 

9.6 

10/7 

7.1 

8.3 

10/8 

6.0 

7.2 

With  the  aid  of  the  above  table,  we  decide  to  use  subarrays  of  8  by  I  elements,  horizontally 
orientated  and  with  an  element  spacing  around  0.73.. 

6.2  Rectangular  grid  with  matching  structure 

The  antenna  geometry  is  shown  in  figure  18;  the  dimensions  obtained  by  optimisation  are  stated 
in  table  4: 


I 
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Figure  18:  Antenna  geometry 

Table  4:  Dimensions  rectangular  grid  antenna  with  matching  structure 


The  worst  case  return  loss  for  the  two  orthogonal  polarisations  as  function  of  the  frequency  over 
the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  19. 

The  worst  case  polarisation  isolation  for  the  two  orthogonal  polarisations  as  function  of  the 
frequency  over  the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  20. 

The  150  MHz  bandwidth  (operational  with  margins)  is  indicated  with  vertical  lines. 
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RETURN  LOSS  (dB) 


FREQUENCY  (GHz) 

—  TE-10  +  TE-01  *  SPECIFICATION 


Figure  19: 


Worst  case  return  loss  over  scan  volume  vs  frequency;  rectangular  grid  with 
matching  structure 
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POLARISATION  ISOLATION  (dB) 


—  TE-10  EXCITED  +  TE-01  EXCITED  *  SPECIFICATION 


Figure  20:  Worst  case  polarisation  isolation  over  scan  volume  vs  frequency;  rectangular 

grid  with  matching  structure 

6.3  Rectangular  grid  without  matching  structure 

The  antenna  geometry  is  shown  in  figure  21;  the  dimensions  obtained  by  optimisation  are  stated 
in  table  5: 


Figure  21:  Antenna  geometry 

Table  5:  Dimensions  rectangular  grid  antenna  without  matching  structure 


S  (mm) 

39.60 

D  (mm) 

35.88 

T  (mm) 

39.60 

T-,  (mm) 

0.00 

A  (mm) 

35.88 

T-,  (mm) 

0.00 

B  (mm) 

35.88 

Er 

C  (mm) 

35.88 

The  worst  case  return  loss  for  the  two  orthogonal  polarisations  as  function  of  the  frequency  over 
the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  22. 

The  worst  case  polarisation  isolation  for  the  two  orthogonal  polarisations  as  function  of  the 
frequency  over  the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  23. 

The  150  MHz  bandwidth  (operational  with  margins)  is  indicated  with  vertical  lines. 
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RETURN  LOSS  (dB) 


♦  TE-10  -t-TE-01  *  SPECIFICATION 


Figure  22:  Worst  case  return  loss  over  scan  volume  vs  frequency;  rectangular  grid  without 

matching  structure 


POLARISATION  ISOLATION  (dB) 


♦TE-10  EXCITED  +TE-01  EXCITED  *  SPECIFICATION 


Figure  23:  Worst  case  polarisation  isolation  over  scan  volume  vs  frequency;  rectangular  grid 
without  matching  structure 
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6.4  Triangular  grid  with  matching  structure 

The  antenna  geometry  is  shown  in  figure  24;  the  dimensions  obtained  by  optimisation  are  stated 
in  table  6: 


Figure  24:  Antenna  geometry’ 

Table  6:  Dimensions  rectangular  grid  antenna  with  matching  structure 


S  (mm) 

37.87 

D  (mm) 

21.47 

T  (mm) 

37.50 

x2  (mm) 

0.00 

A  (mm) 

32.94 

!•,  (mm) 

0.95 

B  (mm) 

£r 

5.58 

C  (mm) 

C2  (degrees) 

69.41 

The  worst  case  return  loss  for  both  orthogonal  polarisations  as  function  of  the  frequency  over 
the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  25. 

The  worst  case  polarisation  isolation  for  both  orthogonal  polarisations  as  function  of  the 
frequency  over  the  scan  volume  (±  1 5°  azimuth,  ±  2°  elevation)  is  shown  in  figure  26. 
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Figure  25:  Worst  case  return  loss  over  scan  volume  vs  frequency;  triangular  grid  with 

matching  structure 


Figure  26: 


Worst  case  polarisation  isolation  over  scan  volume  vs  frequency;  triangular 
grid  with  matching  structure 
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6.5  Triangular  grid  without  matching  structure 

The  antenna  geometry  is  shown  in  figure  27;  the  dimensions  obtained  by  optimisation  are  stated 
in  table  7; 


B 


radiating  aperture 


Legend: 

S,T  Grid 

A,B  Waveguide 

C.D  Iris 


Figure  27:  Antenna  geometry 

Table  7:  Dimensions  rectangular  grid  antenna  with  matching  structure 


S  (mm) 

39.60 

D  (mm) 

37.00 

T  (mm) 

43.84 

t7  (mm) 

0.00 

A  (mm) 

37.00 

T-,  (mm) 

0.00 

B  (mm) 

37.00 

er 

C  (mm) 

37.00 

£2  (degrees) 

62.76 

The  worst  case  return  loss  for  both  orthogonal  polarisations  as  function  of  the  frequency  over 
the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  28. 

The  worst  case  polarisation  isolation  for  both  orthogonal  polarisations  as  function  of  the 
frequency  over  the  scan  volume  (±  15°  azimuth,  ±  2°  elevation)  is  shown  in  figure  29. 


i 


Figure  28:  Worst  case  return  loss  over  scan  volume  vs  frequency;  triangular  grid  without 

matching  structure 


Figure  29: 


Worst  case  polarisation  isolation  over  scan  volume  vs  frequency;  triangular 
grid  without  matching  structure 
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With  the  results  presented  in  chapter  6  as  basis  we  can  design  a  subarray  to  be  used  in  the  SAR 
antenna.  In  order  to  reduce  the  costs  of  the  antenna,  phase  shifting  on  element  to  element  basis  is 
not  recommended.  Due  to  the  limited  scan  volume,  elements  can  be  grouped  in  subarrays. 

In  the  following  the  subarray  design  method,  using  the  software  described  in  the  foregoing  will 
be  described  and  finally  some  subarray  designs  will  be  presented. 

From  now  on,  triangular  lattices  will  not  be  taken  into  account  any  more  because  of  the 
problems  that  wili  arise  due  to  the  gaps  that  will  appear  at  the  tile  boarders. 

7.1  Subarray  design  strategy 

Figure  30  shows  an  array  composed  of  8x1  subarrays.  Elements  in  the  centre  subarray  are 
numbered  from  1  to  8: 


¥ 

0 


□□□□□□□□ 

□□□□□□□□ 

□  □□:□□□□□ 

□□□□□□□□ 

□□□□□□□□ 

□□□□□□□□ 

□ □□□□□□□ 

□□□□□□□□ 

□  □□□□□□□ 

□  □□!□□□□□ 

□□□□□□□□ 

□□□□□□□□ 

□  □  □  □  □  □ 

□□□□□□00 

□  □□,□□□□□ 

□□□□□□□□ 

□□□□□□□□ 

□□□□□□□□ 

0 

0 

JL 


Figure  30:  Centre  subarray  in  array  environment 

To  calculate  the  centre  subarray  reflection  coefficient,  three  steps  need  to  be  performed: 

•  1.  Calculate  the  coupling  coefficients  for  a  finite  set  of  elements  using  infinite  array  data. 

•  2.  Recalculate  the  reflection  coefficient  for  element  i;  i  =  1,2, ...,  8.  In  other  words:  Sum  the 

coupling  coefficients  within  a  window  around  each  element  with  the  appropriate  subarray 
phase  factors.  (Figure  18  shows  a  7x7  window'  for  element  8). 

•  3.  Combine  the  so  obtained  element  in  subarray  reflection  coefficients  to  get  the  overall 

subarray  reflection  coefficient. 
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Step  1  is  obvious:  Apply  a  two  dimensional  Fourier  transformation  on  the  infinite  array 
reflection  coefficients  to  obtain  a  finite  set  of  coupling  coefficients.  Step  2  needs  some  further 
explanation: 

Since  phaseshifters  are  applied  on  subarray  level,  equation  (3.5.2)  needs  to  be  modified  The 
element  in  subarray  reflection  coefficient  is  now  given  by: 


r,(^o»<Po)=  X  =  . 8 

m=-M  n=-N 

with: 

^,  =  ^,(8  dx).msub 
^  =  %.dynsub 

and: 

T,  =  £0sin(i30)cos((p0) 

%-k0  sin(t30)sin(cp0) 


Further: 

nsub  =  n 


(7.1.1) 


(7.1.2a) 


(7.1.2b) 


(7.1.3a) 


(7.1.3b) 


(7.1.3c) 


(7.1.4a) 


a 


msub  =  \ 


c 


if  (m  >  0)A|  a  <  <  b 

+  n 

T 


if  (w  ^  0)aI  c  <  — — —  <  d 


| ,a,be  N 
\\c,d  €  N 


(7.1 ,4b) 
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With  this  design  strategy,  the  subarray  reflection  coefficient  for  the  element  presented  in 
paragraph  6.2,  using  an  ideal  combining  network,  was  calculated  as  function  of  scan  angles  0O 
and  <p0  for  a  frequency  of  5.3  GHz.  A  63x63  window  was  applied.  The  results  are  shown  in 
figure  3 1 : 


Figure  31:  Amplitude  subarray  reflection  coefficients  vs  scan  angles,  f  =  5.3GHz,  TEI0- 

mode  polarisation 

As  a  check,  comparison  with  figures  8,9,10  shows  that  the  reflection  coefficient  at  broadside  is 
the  same  as  for  an  element  in  an  array  scanned  on  element  to  element  basis.  This  is  in  agreement 
with  the  expectations. 

Also  note  that  the  reflection  coefficient  in  the  plane  <p0=90  degrees  is  the  same  as  the  reflection 
coefficient  shown  in  figure  8  where  phaseshifters  were  applied  on  element  to  element  basis.  This 
too  is  in  agreement  with  the  expectations,  since  in  the  plane  <p0=90  degrees  nothing  has  changed 
due  to  the  width  of  the  subarray  (1  element).  Of  course  the  behaviour  in  the  planes  <p0=0  and  <p 
0=45  has  changed  in  respect  with  the  situation  in  which  phaseshifters  were  applied  on  element  to 
element  basis. 


In  order  to  demonstrate  the  influence  of  the  coupling  coefficients,  subarray  reflection 
coefficients  as  function  of  the  scan  angles,  for  this  particular  design,  are  shown  in  figure  32  for 
different  applied  window  sizes: 
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Figure  32: 


Amplitude  subarray  reflection  coefficients  vs  scan  angles  for  different  applied 
window  sizes;  f  -  5.3GHz,  TEl0-mode  polarisation 
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The  above  figure  shows  again  that  adding  more  elements  leads  to  a  higher  fluctuation  frequency 
and  a  lower  fluctuation  amplitude 

The  figure  also  shows  the  importance  of  taking  enough  coupling  coefficients  in  consideration 
The  figure  gives  the  impression  that  when  only  the  coupling  coefficients  within  the  centre 
subarray  are  taken  into  account,  highly  erroneous  reflection  coefficients  will  result. 

To  test  this  assumption,  the  centre  subarray  reflection  coefficients  a?  function  of  the  scanangles 
are  calculated  taking  into  account  only  the  coupling  coefficients  within  this  subarray.  The 
results,  together  with  the  results  when  a  63x63  window  is  applied  are  shown  in  figure  33. 

The  results  clearly  show  the  validity  of  ti.2  above  stated  assumption.  Thus. 

One  can  not  design  a  phased  array  antenna  that  is  scanned  on  subarray  level,  taking  into 

account  only  the  isolated  subarrav  behaviour,  even  when  this  isolated  subanax  shows  the 

same  behaviour  as  one  isolated  element.  The  element  to  element  coupling  mav  be  low .  but  due 

to  the  size  ot  the  complete  phased  array  antenna,  inclusion  of  total  array  coupling  behaviour 

is  of  paramount  importance ! 

To  demonstrate  that  the  above  statement  is  more  or  less  radiator  type  independent,  we 
performed  a  similar  analysis  on  a  microstrip  patch  array  antenna. 

With  coupling  coefficients  of  a  7x7  single  polarisation,  single  layer,  probe  fed  microstrip  patch 
array  antenna,  we  calculated  the  reflection  coefficient  of  an  embedded  4x1  subarray  and  the 
reflection  coefficient  of  an  isolated  4x  1  subarray''.  The  results  are  shown  in  figure  34  for  a 
frequency  of  3  GHz.  Element  distances  are  50  mm,  the  thickness  of  the  dielectric  substrate 
(er=2.94)  is  6  mm. 

Comparison  of  figures  33  and  34  shows  that  for  an  element  optimised  in  an  infinite  array 
environment,  the  embedded  subarray  reflection  coefficient  is  better  than  the  isolated  subarray 
reflection  coefficient,  proving  the  validity  of  our  design  strategy. 


1  These  coupling  coefficients  were  kindly  supplied  by  A.B.  Smolders.  Electromagnetics  Division, 
Faculty  of  Electrical  Engineering,  Eindhoven  University  of  Technology,  P.O.  Box  513,  5600  MB 
Eindhoven,  The  Netherlands. 


TNO-report 


Page  FEL-94-B062 

68 


Figure  33:  Amplitude  waveguide  subarray  reflection  coefficients  vs  scan  angles,  with  and 

without  taking  into  account  the  array  environment;  f  =  5.3GHz,  TEl0-mode 
polarisation,  8x1  subarray 


I 

I 


TNO-reporl 


Figure  34: 


Amplitude  patch  antenna  subarray  reflection  coefficients  vs  scan  angles,  with 
and  without  taking  into  account  the  array  environment;  f  =  3GHz,  TEj0-mode 
polarisation,  4x1  subarray 
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7.2  Subarray  designs 

Because  in  the  subarray  designs  to  be  presented,  co  and  cross  polarisation  radiation  patterns  will 
be  shown,  we  start  with  a  brief  discussion  on  the  calculation  of  co  and  cross  polarisation 
behaviour.  These  radiation  patterns  are  based  on  radiation  patterns  of  a  singly  excited  element  in 
an  infinite  array  environment.  We  present  this  discussion,  because  of  the  fact  that  calculating  the 
radiation  pattern  of  a  singly  excited  element  in  an  infinite  array  environment  -  bearing  in  mind 
that  our  analysis  program  is  based  on  a  uniform  excited  infinite  array  -  is  not  al  all  that  obvious. 

7.2.1  Co  and  cross  polarisation  analysis  method 

Use  will  be  made  of  the  fact  that  by  taking  a  Fourier  series  transform  of  the  field  solution  for  the 
infinitely  excited  array  one  may  obtain,  uniquely,  the  field  solution  for  the  singly  excited 
element  in  the  environment  of  the  array  [26].  More  precisely:  If  a  Fourier  analysis  is  made  of  the 
aperture  fields  of  an  active  infinite  array,  which  by  Floquet's  theorem  are  periodic  functions  of 
the  steering  angles,  the  mth  Fourier  coefficient  is  then  the  aperture  fields  across  the  mth  element 
when  the  zerorh  element  is  excited  alone  [27). 

7.2.1.1  Aperture  fields 

Started  is  with  the  Huygen's  principle,  see  figure  35: 


(Electric  and  Magnetic  fields  desired) 


Actual  sources 

0- 


Figure  35:  Radiation  through  an  aperture;  Huygens'  principle 

The  tangential  components  of  the  electric  and  magnetic  fields  in  the  aperture  need  to  be  known. 
Although  these  fields  arise  from  sources  to  the  left,  they  may  be  considered  to  be  produced  by 
equivalent  sources  located  in  the  aperture  plane.  The  fields  are  described  by  [28]; 
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£  =  -;(o4-^V(V.a)--Vx£  (7.2.1) 

k  E 

with: 

e~>k'  _ 

A  =  H~ — N  (7.2.2a) 

4rtr 

e~ikr  - 

£  =  e- — £  (7.2.2b) 

47tr 


=  f  Jse+J*rcos,r^5  (7.2.3a) 

Z  =  j  Mxe+jkr'™lT)dS  (7.2.3b) 

5 


J,  =nxH 


(7.2.4a) 


A/(  =  -«:r£ 


(7.2.4b) 


with  £  the  electric  field,  A  and  £  retarded  magnetic  and  electric  vector  potentials,  L  and  N 

magnetic  and  electric  radiation  vectors,  Ms  and  Js  surface  magnetic  and  electric  currents  and  k 
the  wavenumber: 


(7.2.5) 


n  is  a  unit  vector  normal  to  the  aperture  plane,  r'  and  4*’  are  defined  in  figure  36  [28]: 
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Figure  36:  Coordinates  of  a  general  current  element  in  the  volume  V  and  a  distant  point 

Q.  where  fields  are  to  be  evaluated,  with  respect  to  an  origin  of  coordinates 
near  the  sources 

7.2.1.2  Infinite  waveguide  array 

When  the  waveguide  properties  are  described  in  a  Cartesian  x,y,z-coordinate  system,  with  the 
wave  propagation  directed  along  the  z-axis,  equation  (7.3.4a)  can  be  written  as: 


Js  =  hxH  =  -Hxet  +  Hxev 


(7.2.6) 


The  tangential  magnetic  field  strengths  (the  Fourier  coefficients  mentioned  in  paragraph  7.2.1) 
are  given  by  equations  (3.3.27)  and  (3.3.28): 


=(-y!4,W«p07'(4,+/v  +we,e0vwr(4)+M)xw  lkTJ  (7.2.7a) 

(7.2.7b) 


with: 
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=  e 


(7.2.8) 


and  the  other  quantities  defined  in  chapter  3  and  [12]. 

Equation  (7.2.8)  is  a  phased  array  term,  related  to  the  element  positions  in  the  array,  and  is  only 
valid  for  discrete  x  and  y. 

Now,  we  want  to  calculate  the  radiation  pattern  of  the  singly  excited  unit  cell  in  the  infinite 
array  environment.  Since  antennas  are  designed  such  that  only  one  grating  mode  is  radiating  in 
free  space,  we  take  p=q=0.  Since  we  look  at  the  unit  cell,  we  take  x=y=0.  This  yields. 


Hx= 1 

(— Y(4)<»w007'l4'*oo  +  mr4eoVooT^oo)  /  *,n2 

(7.2.9a) 

H,  =  1 

(-Y,4,oo  VooT- ,,+oo  -  mre0umf(4Hoo )  / 

(7.2.9b) 

Since  our  origin  is  placed  in  the  aperture  plane,  we  take  r'=4/,=0  in  equation  (7.2.3a).  With  unit 
cell  area  s.f.sin(£2)  (see  chapter  3)  and  use  of  equations  (7.2.2a),  (7.2.3a),  (7.2.6)  and  (7.2.9): 


e~lkr 

A  -  M(,.s./.Sin(£2)— — 
47tr 


Y  «>v(io*  +  coe  e0«()0T,4K,«»  „  Y,4,<K.u(K)r4'+.«l-WE,<E(,v(l()r 

- - — - - e - - — 


V 2 


(7.2.10) 


Since  VxF  =  0  and  V  •  A  -  0 ,  equation  (7.3. 1 )  reduces  to: 


E  =  -  70X4 


(7.2.11) 


So,  finally: 


Ex  =-y'ox4,  =-/co[t0.j.t.sin(Q). 


e  Jkr  y^ooVooT  oo  +  m^u^r^oo 


4nr 


kT 

‘(t) 


(7.2.12a) 


£,  =  -  yox4v  =  +  y'(D|i0 .  s.  t.  sin(Q). 


e->k'  7M)oomqo7',4H(x>  -CQ£r<e0voor<4Hoo 


4Jtr 


kT 

*nr> 


2 


(7.2.12b) 
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£  =  -  jtiiA.  =  0 


(7.2.12c) 


7.2.1.3  Co  and  cross  polarisation 

Applying  Ludwig's  third  definition  [29]  and  omitting  the  terms  exp(-jkr)/4jtr,  we  get,  after 
normalisation: 


Eu,Voi  =+(l-cos(T>))sin((p)cos((p).-^7.(y'4'ooi’o07''4H«.  +  a)eri£oW0o7''4'+u  <- 

MKi 

{l  -  sin2  (cp)(l  -  cos($))}.— y  .(y<4*ooh007'|4)+oo  -  coe,4E0v007r(4)"oo ) 

ki„ 

(7.2.13a) 


Krp.,1  -  -{l-COS'(<p)(l-COs(d) — -.(y^Woo^'Vi  +COEr4£0W007''',KOo)  + 

^7■,<, 

-(l-cos(v))sin(cp)cos(9).-^T.(yl4,ooM007’'4Hoo-o)£r4E0i'oorJ|*oo) 

kTm 


(7.2.13b) 

Figure  37  shows  measured  and  calculated  active  element  co  and  crosspolarisation  radiation 
patterns  for  the  EXPAR  test  antenna  (see  chapter  4)  in  the  worst  case  plane  (f  =  f2,  plane  (p=55 
degrees): 


t  I 
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Amplitude  (dB) 


Theta  (degrees) 
*  calc.  meas. 


Figure  37:  EXPAR  worst  case  plane  active  element  co  and  cross  polarisation  radiation 

patterns 

Note  the  fluctuations  in  the  measured  patterns  due  to  the  finiteness  of  the  array  (255  elements), 
as  discussed  before. 

7.2.2  Subarray  with  matching  structure 

For  a  8x  1  subarray  with  dimensions  as  stated  in  paragraph  6.2,  the  worst  case  return  loss  over 
the  scanvolume  as  function  of  frequency  is  calculated,  following  the  strategy  explained  in 
paragraph  7.2.  The  results  for  the  two  orthogonal  modes  are  shown  in  figure  38: 
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RETURN  LOSS  (dB)  < 


TE-10  +TE-01  SPECIFICATION 


Figure  38:  Subarray  worst  case  return  loss  over  scan  volume  vs  frequency;  rectangular 

grid  with  matching  structure 

Due  to  the  very  low  values  of  the  coupling  coefficients  for  the  not  excited  mode,  it  is  not 
possible  to  calculate  the  subarray  polarisation  isolation  in  the  same  way  as  outlined  in  paragraph 
7.2.  Because  of  the  agreement  between  worst  case  return  loss  for  the  subarray  (above  figure)  and 
worst  case  return  loss  for  the  element  in  an  infinite  array  environment  (see  paragraph  6.2,  figure 
19),  it  is  allowed  to  conclude  that  the  polarisation  isolation  will  be  comparable  to  that  of  the 
element  in  an  infinite  array  environment  (see  paragraph  6.2,  figure  20). 

For  the  same  reason,  it  is  allowed  to  construct  subarray  radiation  patterns  based  on  the  radiation 
patterns  of  an  element  in  an  infinite  array  environment.  Figure  39  shows  the  element  patterns 
and  figure  40  shows  the  subarray  radiation  patterns  based  on  these  element  patterns. 

The  aperture  efficiency  can  be  approximated  in  first  order  approximation  by  [30]: 
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.  Embedded  subarray  gain  at  input  port 

Aperture  efficiency  = - — =-r — — — — - - — - — 

r  4ti.  Subarray  aperture 

4  n 


J.8.S.  T.  cos(  ).  ( 1  -  |r(  d0 ,  tp0  )|^  ) 


471.8.5.7' 


li3„=po=0 


(7.2.14) 


=  Hr(o.<C 


with  S  and  T  the  element  spacings  in  x  and  y  directions,  respectively. 

Ohmic  losses  in  the  subarray  feeding  network  are  not  taken  into  account  in  the  above  equation. 

The  aperture  efficiency  is  shown  in  figure  41.  The  bandwidth  (operational  with  margins)  is 
indicated  with  vertical  lines. 
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Figure  40:  Radiation  patterns  subarray  in  array  environment;  rectangular  grid  with 

matching  structure,  f  =  5.3GHz 
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Figure  41:  Aperture  efficiency  subarray;  rectangular  grid  with  matching  structure,  ideal 

lossless  feeding  network 

7.2.3  Subarray  without  matching  structure 

For  a  8x1  subarray  with  dimensions  as  stated  in  paragraph  6.3,  the  worst  case  return  loss  over 
the  scanvolume  as  function  of  frequency  is  calculated,  following  the  strategy  explained  in 
paragraph  7.2.  The  results  for  the  two  orthogonal  modes  are  shown  in  Figure  42: 
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Figure  42:  Subarray  worst  case  return  loss  over  scan  volume  vs  frequency;  rectangular 

grid  without  matching  structure 

For  this  configuration,  it  also  was  impossible  to  calculate  the  polarisation  isolation,  following 
the  rules  of  paragraph  7.2,  due  to  the  very  low  values  of  the  coupling  coefficients  for  the  not 
excited  mode. 

Because  of  the  agreement  between  worst  case  return  loss  for  the  subarray  (above  Figure)  and 
worst  case  return  loss  for  the  element  in  an  infinite  array  environment  (see  paragraph  6.3,  figure 
22),  it  is  allowed  to  conclude  that  the  polarisation  isolation  will  be  comparable  to  that  of  the 
element  in  an  infinite  array  environment  (see  paragraph  6.3,  figure  23). 

For  the  same  reason,  it  is  allowed  to  construct  subarray  radiation  patterns  based  on  the  radiation 
patterns  of  an  element  in  an  infinite  array  environment.  Figure  43  shows  the  element  patterns 
and  figure  44  shows  the  subarray  radiation  patterns  based  on  these  element  patterns. 

The  subarray  aperture  efficiency  is  shown  in  figure  45. 
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Figure  43: 


Element  radiation  patterns  in  infinite  array  environment;  f  =  5.3GHz, 
rectangular  grid  without  matching  structure,  TEl0-mode  polarisation 
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Figure  44:  Radiation  patterns  subarray  in  array  environment;  rectangular  grid  without 

matching  structure,  f  -  5.3GHz 
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Figure  45:  Aperture  efficiency  subarray;  rectangular  grid  without  matching  structure, 

ideal  lossless  feeding  network 


7.3  Choice  subarray  design 

In  an  earlier  stage  we  already  rejected  designs  with  a  triangular  lattice.  At  this  stage  we  have  to 
make  a  choice  between  the  design  of  paragraph  6.2  (7.2.2)  and  the  design  of  paragraph  6.3 
(7.2.3). 


Although  the  characteristics  of  the  design  pmtnts<Lin  jtuagtmhJiJ.  £7,2,2)  are  Inter,  w? 

:hoose  for  the  design  of  paragraph  6.3  (7.2.3). 


The  reason  for  this  choice  is  that  this  design  is  much  easier  to  construct  (no  irises  and  dielectric 
matching  sheet).  Furthermore,  we  do  not  have  to  deal  with  tolerance  problems  in  applying  the 
dielectric  sheet  and  the  characteristics  of  this  design  still  meet  the  specifications. 
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8  WAVEGUIDE  EXCITER  &  SUBARRAY  FEEDING  NETWORK 


In  the  foregoing  chapters,  we  assumed  that  the  TEjQ-mode  and/or  the  TEQj-mode  were  excited 
in  the  waveguide-  without  specifying  how  this  was  accomplished.  In  calculating  the  subarray 
reflection  coefficient  we  assumed  perfectly  excited  waveguide  modes  and  an  ideal  1  to  8 
combining  network. 

This  chapter  deals  with  the  practical  implementation  of  exciting  the  waveguide  modes  and 
combining  the  elements  in  the  subarray. 

The  waveguide  exciter,  to  be  discussed  in  the  next  paragraph,  is  ment  for  breadboard  purposes 
only.  The  exciter  is  not  optimal  for  space  qualification,  but  a  first  step  in  the  development  of  a 
completely  space  qualified  design.  The  space  qualified  design  will  be  generated  under 
workpackage  2000  [1].  Since  workpackage  1000  only  asks  a  preliminary  space  borne  SAR 
antenna  design,  the  exciter  to  be  discussed  will  do  for  the  moment. 


8.1  Preliminary  waveguide  exciter  design 

In  order  to  make  an  exciter  as  flat  as  possible  -  radiator  thickness  is  an  important  design 
parameter  -  we  use  a  microstrip  patch  structure  to  excite  the  TEjo-mode  and  TEgi-mode. 

This  may  seem  odd  at  first  sight,  since  one  of  the  antenna  options  is  a  microstrip  patch  antenna 
11],  but  -  and  this  can  not  be  emphasised  enough  -  the  patch  structure  in  itself  is  not  used  as  an 
antenna  element. 

We  do  not  have  to  deal  with  microstrip  patch  mutual  coupling,  ground  waves  and  other  effects 
of  microstrip  patch  array  antennas.  The  patch  is  used  as  exciter  of  the  waveguides,  and  the 
behaviour  of  waveguide  antennas  can  be  simulated  with  great  accuracy,  as  is  thoroughly 
demonstrated  in  the  foregoing  chapters. 

The  idea  of  using  a  microstrip  patch  antenna  is  not  new;  results  have  recently  been  published 
[31]. 

Because  of  the  involvement  of  TNO-FEL  in  the  design  and  realisation  of  a  C-band  airborne 
SAR  antenna  (PHARUS:  PHased  ARray  Universal  SAR)  [32],  single  layered  dual  polarised 
probe-fed  patch  elements  were  available  and  these  were  used  as  basis  for  the  waveguide  exciter. 
The  exciter  is  shown  in  figure  46.  For  the  sake  of  clarity,  quarter  lambda  airstripline 
transformers  in  the  coaxial  connectors,  are  not  drawn. 
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Figure  46:  Top  and  side  view  waveguide  exciter 

A  demonstrator  has  been  realised.  The  dimensions,  with  reference  to  the  above  figure,  are  stated 
in  table  8.  Results  of  a  electromagnetic  3D  simulator  as  well  as  measurement  results  are  shown 
in  figure  47a  and  47b.  Measured  and  simulated  are  the  behaviour  of  the  exciter  placed  in  the 
waveguide 


Table  8:  Dimensions  waveguide  exciter  demonstrator 


Ww  (mm) 

34.38 

hx  (mm) 

2.50 

Lw  (mm) 

34.38 

hy  (mm) 

2.50 

Wp  (mm) 

17.50 

h  (mm) 

0.035 

Lp  (mm) 

17.50 

d  (mm) 

1.57 

wx  (mm) 

0.8 

Er 

2.33 

RETURN  LOSS  (dB) 


-*•  MEASUREMENT  +3D  SIMULATION 


Figure  47a:  Measured  and  simulated  waveguide  exciter  return  loss 
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Figure  47b:  Measured  and  simulated  waveguide  exciter  polarisation  isolation 

The  mismatch  between  simulations  and  measurements  in  the  second  half  of  the  frequency  band 
considered  is  most  probably  due  to  differences  in  the  stub  dimensions  (see  figure  25)  in 
simulations  and  measurements.  During  the  measurements,  the  stubs  were  lost  and  replaced  with 
wires  of  length  equal  to  the  stublengths  as  shown  in  figure  46. 

Nevertheless,  the  figures  clearly  show  that  we  are  able  to  simulate  the  behaviour  of  the  exciter 
placed  in  the  waveguide. 

It  is  interesting  to  note  that  the  behaviour  measured  at  the  right  side  of  the  resonance  frequency 
also  shows  up  in  the  simulations.  This  behaviour  must  be  due  to  a  mismatch  of  the  quarter 
lambda  transformers  mentioned  before. 

Furthermore,  the  exciter  satisfies  the  specifications  within  the  150  MHz  frequency  band. 

8.2  Space  qualified  waveguide  exciter  design 

The  foregoing  paragraph  has  shown  that  we  have  the  tools  necessary  to  design  a  microstrip 
patch  waveguide  exciter. 

At  the  moment,  running  ahead  of  workpackage  2000,  we  are  working  on  a  aperture  coupled 
microstrip  patch  waveguide  exciter  that  has  no  soldering  connections  and  will  therefore  be  space 
qualified. 
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8.3  Subarray  feeding  network  design 

In  order  to  demonstrate  our  ability  to  design  a  subarray  feeding  network,  running  ahead  of 
workpackage  2000  and  the  space  qualified  exciter  design,  a  feeding  network  in  microstrip  is 
designed.  Material  parameters  are  given  in  table  9: 


Table  9:  Subarray  feeding  network  material  parameters 


|  PARAMETER 


VALLE 


fr 

2.33 

tan(5) 

0.0012 

substrate  height 

0.4  mm 

me'al  thickness 

0.035  mm 

metal  rho 

1 .25 

metal  roughness 

0 

Figures  48a,  b  and  c  show  the  patches  in  the  waveguides,  the  feeding  network  and  the 
combination  of  these,  respectively: 


Figure  48a  Subarray:  exciters  withtn  waveguides 


Figure  48b:  Subarray:  microstrip  feeding  network 
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Figure  48c:  Subarray:  exciters  within  waveguides  and  microstrip  feeding  network 

Note  that  the  feeding  network  completely  fits  in  the  subarray. 

The  behaviour  of  this  feeding  network  has  been  simulated.  The  results  are  shown  in  figures  49, 
51  and  52. 

These  figures  show  that  all  specifications  are  met, except  for  the  loss  goal  of  ldB/m. 

With  a  subarray  length  of  31.68  cm  (design  of  paragraph  6.2),  figure  49  indicates  that  the  loss 
will  be  1.19dB/m.  Since  this  feeding  network  is  a  preliminary  design,  we  are  convinced  that  we 
can  design  a  feeding  network  that  even  meets  this  rather  tight  specification.  This  can  be 
accomplished  by  using  a  thicker  substrate.  A  thickness  of  1  mm  seems  to  be  possible  without 
getting  problems  like  surface  waves. 


Shown  is  that  we  have  the  tools  and  experience  to  design  a  subarray  feeding  network. 
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EEsof  -  liDra  -  wed  Jan  12  09:48:03  1994  -  second 
c  t)B[  S21 ! 


Figure  51:  Transmission  coefficient  1:8  subarray  power  divider 
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9  RADIATING  PLANE  THICKNESS 


In  order  to  estimate  the  depth  of  the  waveguides,  we  have  to  look  at  the  non-propagating 
waveguide  modes.  These  modes  are  introduced  at  the  exciter  side  of  the  waveguide  and  at  the 
aperture  (iris)  side  of  the  waveguide.  The  length  of  the  waveguide  must  be  chosen  such,  that 
these  non  propagating  higher  order  modes  -  the  waveguides  can  only  propagate  the  TEjQ-mode 
and  the  TEQj-mode  -  are  sufficiently  attenuated,  so  that  no  higher  order  mode  interaction 
between  exciter  side  and  aperture  side  can  occur. 


9.1  Attenuation  constant 


Since  waves  are  propagating  or  decaying  with  a  factor  exp(— yy1  Umnz),  z  being  the  distance  in 
the  waveguide  direction,  the  attenuation  constant  for  non  propagating  modes  is  given  by 

a  =  -Real(-yy<l)mn)- 

The  attenuation  constants  for  the  waveguide  modes  of  the  design  presented  in  paragraph  6.2  at 
the  centre  frequency  are  stated  in  table  10  and  figure  52.  These  attenuation  constants  change 
very  little  over  the  frequency  band  of  150  MHz. 


Table  10:  Attenuation  constants  higher  order  modes  at  5.5  GH: 


WAVEGUIDE  MODE 


a  (nr1) 


a  (dB/m) 


TE-02 

155.12 

43.81 

TE-11 

76.62 

37.69 

TE-12 

182.09 

45.21 

TE-20 

155.12 

43.81 

TE-21 

182.09 

45.21 

TE-22 

245,86 

47.81 

TE-30 

263.71 

48.42 

TE-31 

280.43 

48.96 

TE-32 

325.46 

50.25 

TE-40 

364.99 

51.25 

TE-41 

377.24 

51.53 

TE-42 

411.82 

52.29 

TE-50 

463.77 

53.33 

TE-51 

473.47 

53.51 

TE-52 

501.46 

54.00 

TM-11 

76.62 

37.69 

TM-12 

182.09 

45.21 

TM-21 

182.09 

45.21 

TM-22 

245.86 

47.81 

TM-31 

280.43 

48.96 

TM-32 

325.46 

50.25 

TM-41 

377.24 

51.53 

TM-42 

411.82 

52.29 

TM-51 

473.47 

53.51 

TM-52 

501.46 

54.00 
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Figure  52:  Higher  order  mode  attenuation  constants  at  5.3GHz 

Table  6  and  figure  52  show  that  the  TE]j-mode  and  TM||-mode  are  the  modes  that  will 
determine  the  waveguide  hight. 

9.2  Attenuation  levels 

Attenuation  levels  as  function  of  the  waveguide  height  (  exp(az)  )  are  given  in  table  1 1  and 
figure  53. 
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Table  II:  Attenuation  levels  ( dB )  as  function  of  waveguide  thickness 

MODE  z=lcm  z  =  2cm  I  z  =  3cm  z  =  4cm 


TE-02 

13.48 

26.95 

40.43 

53.90 

TE-11 

6.66 

13.31 

19.97 

26.62 

TE-12 

15.82 

31.64 

47.45 

63.26 

TE-20 

13.47 

26.94 

40.42 

53.89 

TE-21 

15.82 

31.63 

47.45 

63.27 

TE-22 

21.35 

42.71 

64.06 

85.42 

TE-30 

22.91 

45.81 

68.72 

91.62 

TE-31 

24.36 

48.72 

73.08 

97.43 

TE-32 

28.27 

56.54 

84.81 

113.08 

TE-40 

31.70 

66.28 

95.11 

126.81 

TE-41 

32.76 

65.53 

98.30 

131.06 

TE-42 

35.77 

71.54 

107.31 

143.08 

TE-50 

40.31 

80.57 

120.85 

161.13 

TE-51 

41.12 

82.25 

123.37 

164.50 

TE-52 

43.55 

87.11 

130.66 

174.22 

TM-11 

6.66 

13.31 

19.97 

26.62 

TM-12 

15.82 

31.64 

47.45 

63.26 

TM-21 

15.82 

31.64 

47.45 

63.26 

TM-22 

21.35 

42.71 

64.06 

85.42 

TM-31 

24.36 

48.72 

73.08 

97.43 

TM32 

28.27 

56.54 

84.81 

113.08 

TM-41 

32.76 

65.53 

98.30 

131.06 

TM-42 

35.77 

71.54 

107.31 

143.08 

TM-51 

41.12 

82.25 

123.37 

164.50 

TM-52 

43.55 

87.11 

130.66 

174.22 
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Figure  53:  Attenuation  levels  vs  radiator  thickness  atf  =  5.3GHz 


9.3  Depth  waveguide  radiators 

Taking  into  account  that  the  higher  order  mode  amplitudes  will  be  small  compared  to  the  excited 
TE]0-mode  and  TEQj-mode  amplitude  coefficients  and  that  we  have  neglected  attenuation 
arising  from  the  finite  conductivity  of  the  waveguide  walls,  looking  at  the  TEjj-mode  and 
TM]  | -mode,  we  choose  a  radiator  thickness  of  2.5  cm.  leaving  0.5  cm  for  the  subarray  feeding 
network./. 

Because  in  all  our  antenna  designs  (see  chapter  6  and  7)  the  waveguide  dimensions  are  more  or 
less  the  same,  it  is  permitted  to  maintain  this  waveguide  depth  for  all  designs. 

Since  the  waveguides  will  be  self-supporting,  the  support  stucture,  common  in  printed  antenna 
SAR  antennas,  is  cancelled,  see  figure  54  [33]  and  figure  55: 


/This  in  agreement  with  the  3cm  for  radiator,  feed  network  and  support  structure  settled  at  the  progress 
meeting  held  at  Dornier  GmbH,  Friedrichshafen  Germany,  November  18  1993. 
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Figure  54: 


Cross  section  of  a  patch  antenna  tile  and  waveguide  antenna  tile 
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Figure  55: 


3D  impression  part  of  waveguide  antenna 
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10  CONCLUDING  REMARKS 


In  the  foregoing  it  has  been  shown  thoroughly  that  at  TNO-FEL  we  have  all  the  tools  necessary 
to  design  a  SAR  waveguide  phased  array  antenna: 

•  The  antenna  can  be  analysed  and  synthesised 

•  Finiteness  effects  can  be  taken  into  account 

•  The  subarray  in  the  array  environment  can  be  analysed;  it  is  demonstrated  that  it  is  of  vital 
importance  to  include  the  array  environment  in  the  subarray  analysis 

•  A  microstrip  patch  waveguide  exciter  can  be  designed  and  analysed 

•  A  microstrip  subarray  feeding  network  can  be  designed  and  analysed 

Different  waveguide  SAR  antenna  designs  are  presented  which  all  satisfy  the  specifications. 

All  the  software  tools  are  available  to  analyse  the  complete  subarray  from  feeding  network  to 
exciter  to  waveguides  all  together. 


In  the  near-future  we  will  perform  measurements  on  a  dual  polarised  (square)  waveguide  array 
antenna.  This  antenna,  shown  in  figure  56,  is  now  under  construction.  _ 
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Figure  56: 


Dual  polarised  waveguide  array  test  antenna 
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APPENDIX  A:  HELPFUNCTIONS  USED  IN  MATCHING  EQUATIONS 


In  the  determination  of  the  helpfunctions  Fj  to  Fjg,  use  is  made  of: 


NE(n)  =  2.sin(— )/ n 
2 


(A. la) 


1,  m  =  n  =  0 
NC(m,  «)  =  •—,  m  =  n*0 
0,  m?sn 


(A. lb) 


0,  m  =  n  =  0 
NS(m,n)  =  -  m  =  n  *  0 
0,  m  *  n 


(A.lc) 


FCT(m,n)  =  njJm  -  12 

Z 


(A. Id) 


FST{m,n)  =  In 


(A.le) 


CCT(m,n,  t|)  = 


7icoslfj4(T)/2;m=a2’4’- 

=  1,3,5,... 


(A. If) 
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SST(m,n,f\ ) 


=  0,2,4,... 
=  1,3,5,... 


With  these  Bessel  function  expanded  integrals: 


(A.lg) 


(A. 2a) 


t.sin(Q)N£([v/1(/ 


F}  =  a.  NC(m,m  ) 

F,  =  b.NS(n,ri  ) 

=  a.NS(m,m  ) 

Fb  =  b.NC{n,ri ) 

F1  =  c.FCT(l,up0c) 

F8  =  FCT{w,vp.q.d ) 


(A. 2b) 

(A.2c) 

(A. 2d) 

(A.2e) 

(A.2f) 

(A.2g) 

(A.2h) 


F9=c.FST{l,up.0c) 


(A.2i) 
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